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Results  are  presented  with  analysis  and  discussion  related 
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FOREWORD 


This  report  was  prepared  by  the  Materials  and  Structures 
Laboratory,  Manufacturing  Development,  General  Motors  Corpora¬ 
tion,  General  Motors  Technical  Center,  Warren,  Michigan,  under 
Contract  DAAG46-69-C-0127  from  the  Army  Materials  and  Mechanics 
Research  Center,  Watertown,  Massachusetts.  The  work  was  per¬ 
formed  during  the  period  30  June  1969  to  15  May  1970  and  was 
technically  monitored  and  administered  by  Mr.  J.  Dignam  and 
Dr.  S.  C.  Chou  of  AMMRC . 


The  objective  of  this  program  was  to  examine  the  behavior 
of  three  aluminum  alloys,  2024-T3,  2014-T6,  and  6061-T6,  in 
tests  related  to  those  needed  to  evaluate  the  vulnerability 
of  interceptor  missiles  by  carrying  out  the  following  programs: 


TASK  A  Study  the  elastic  behavior  of  2014-T6 
aluminum  using  ultrasonic  techniques 
for  temperatures  up  to  400°F.  The 
effect  of  pressure  up  to  8  kilobars 
was  investigated. 

TASK  B  Study  the  yield  surface  of  2024-T3 

aluminum  by  biaxially  stressing  thin 
wall  tubular  specimens .  The  effect 
of  prestraining  on  the  subsequent 
yield  surface  was  investigated  in 
conjunction  with  a  contract  from 
the  Sandia  Laboratories,  No.  53-0204. 

TASK  C  Study  the  high  heating  rate  properties 
of  2014-T6  and  6061-T6  aluminum  at 
temperatures  up  to  700°F.  Uniaxial 
stress  tensile  tests  were  conducted 
after  heat  up  at  rates  up  to 
3  X  10^°F/second. 

TASK  D  Hugoniot  equations  of  state  and  wave 
profile  tests  were  conducted  at  room 
temperature  on  2014-T6  aluminum. 


TASK  E  Spallation  threshold  measurements  on 
2014-T6  aluminum  were  conducted  at 
room  temperature. 


This  report  is  divided  into  four  sections  based  on  the 
type  of  information  discussed.  Section  I  discusses  the  bi¬ 
axial  subsequent  yield  study.  Section  II  discusses  the  high 
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heating  rate  results,  and  Section  III  discusses  the  shock  wave 
data  including  elastic  constants,  hugoniot  and  wave  profile 
results,  and  spallation  results.  Section  IV  contains  conclu¬ 
sions  drawn  from  the  program. 

The  program  was  supervised  and  managed  by  Mr.  S.  J.  Green. 
Principal  investigators  were  Mr.  S.  G.  Babcock  and  Mr.  W.  M. 
Isbell  while  project  scientists  included  Messrs.  J.  D.  Leasia, 
J.  J.  Langan,  F.  L.  Schierloh,  D.  B.  Norvey,  D.  R.  Christman 
and  T.  E.  Michaels. 

The  authors  would  like  to  acknowledge  the  following  in¬ 
dividuals  for  their  competent  assistance  in  the  fulfillment 
of  this  program:  Mr.  M.  C.  Klewicki  and  Mr.  L.  A.  Seitz  for 
their  work  on  the  biaxial  testing,  Mr.  J.  E.  Bonner  for  his 
work  cn  the  high  heating  rate  testing,  and  Mr.  J.  R.  Havens 
and  Mr.  C.  E.  Woodcock  for  their  work  in  the  sho.k  wave 
testing.  In  addition,  special  recognition  is  ex--  nded  to 
Dr.  Arfon  K.  Jones  for  his  technical  comments  and  assistance 
in  this  program. 
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INTRODUCTION 


In  the  vulnerability  assessment  of  reentry  and  interceptor 
vehicles,  the  response  of  the  vehicle  materials  and  structure 
to  stress  wave  and  momentum  loading  must  be  determined.  An 
analysis  of  stress  wave  propagation  is  carried  out  first  to 
determine  the  kinetic  energy  distribution  among  the  fractured, 
spalled,  and/or  debonded  materials  in  the  structure.  Following 
this  is  an  analysis  of  loading  on  the  structure  due  to  the 
momentum  imparted  by  the  stress  wave  interaction.  Computer 
codes  are  used  to  assist  in  the  mathematical  simulation  of 
these  responses,  and  laboratory  or  underground  tests  are  em¬ 
ployed  to  check  the  validity  of  the  computations. 

Inputs  to  these  codes  require  material  properties  in  the  form 
of  constitutive  equations  derived  from  controlled  laboratory 
experiments.  This  program  was  directed  towards  determining 
the  material  properties  required  in  the  computer  codes  men¬ 
tioned  above.  Four  areas  of  experimental  work  were  covered. 

Three  areas,  elastic  property  measurements  using  ultrasonic 
techniques,  shock  wave  experiments  to  determine  hugoniot  equa¬ 
tion  of  state,  complete  wave  profiles  and  spallation  thresholds, 
as  well  as  subsequent  yield  surface  determination  are  familiar 
to  investigators  and  have  been  used  in  past  studies.  One  area, 
high  heating  rate  response,  was  not  well  developed  prior  to 
the  beginning  of  this  program,  and  techniques  had  to  be  developed 
to  the  stage  of  obtaining  valid  data  on  the  materials  of  interest. 

Section  I  deals  with  the  determination  of  the  virgin  and  subse¬ 
quent  yield  surfaces  of  2024-T3  aluminum,  as  determined  by  quasi¬ 
static  biaxial  stress  tests  yielding  data  in  the  four  quadrants 
of  the  axial  stress-circumferential  stress  plane.  The  General 
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Motors  Biaxial  Strain-Rate  Machine  is  described,  along  with 
the  corresponding  variable  input  and  data  reduction  operations. 
Specimen  configurations  are  detailed,  and  the  prestraining  pro¬ 
cess  for  the  subsequent  yield  specimens  is  outlined.  Following 
the  data  presentation  is  a  short  discussion  of  the  literature 
and  a  comparison,  which  leads  to  the  conclusions  for  the  section. 

Section  II  describes  the  technique  developed  to  study  the  heat¬ 
ing  rate  and/or  time-at-temperature  behavior  of  two  aluminum 
alloys:  6061-T6  and  2014-T6.  Temperatures  up  to  700°F  were 

attained  using  a  direct  resistance  heating  technique  at  heat- 

3 

ing  rates  up  to  2.5  X  10  °F/second.  Constant  strain  rate 
10/'sec)  uniaxial  stress  tensile  tests  were  conducted  to 
determine  yield  and  ultimate  strengths  at  various  test  tempera¬ 
tures  and  heating  rates. 

Section  III  describes  the  results  of  measurements  made  on 
2014-T6  aluminum  to  determine  the  elastic  constants  up  to 
400°F  and  8  kilobars  pressure.  Ultrasonic  techniques  were  em¬ 
ployed  to  gather  this  information.  In  addition,  room  tempera¬ 
ture  hugoniot  equation  of  state  and  several  complete  wave  pro¬ 
file  measurements  were  made  and  the  spallation  threshold  was 
determined. 

Section  IV  contains  the  conclusions  drawn  from  these  areas  of 
study.  An  Appendix  is  provided  to  describe  the  physical  and 
chemical  makeup  of  the  three  aluminum  alloys. 
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SECTION  I 

SUBSEQUENT  YIELD  SURFACE  OF 
2024-T3  ALUMINUM 


INTRODUCTION 

j 

The  effect  of  previous  loading  history  on  a  material's  subse¬ 
quent  behavior  has  been  a  subject  of  great  interest.  In  this 
study  the  General  Motors  Biaxial  Machine,  a  unique  dynamic 
loading  apparatus,  was  employed  in  obtaining  the  virgin  and 
subsequent  yield  surfaces  of  2024-T3  aluminum.  The  subsequent 
yield  data  was  gathered  following  the  application  of  an  8% 
axial  compressive  strain  to  the  material.  A  description  of 
the  machine,  specimens,  and  data  handling  techniques  are  in¬ 
cluded  in  this  section.  The  data  itself  is  reviewed  in  re¬ 
gard  to  the  effect  of  the  prestraining  on  the  behavior  of 
the  yield  surface. 


BIAXIAL  TESTING  TECHNIQUE 


The  General  Motors  Biaxial  Strain-Rate  Machine, ^  pictured 
in  Figure  1  and  shown  schematically  in  Figure  2,  develops  a 
biaxial  state  of  stress  by  independently  applying  an  axial 
load  and  a  circumferential  load  (by  fluid  pressure}  on  a 
tubular  specimen  through  two  independent  gas -operated  cylin¬ 
ders.  Strain  rates  from  10-4  to  102/sec  are  obtainable  and 
may  be  held  constant  within  a  factor  of  two  for  the  elastic 
range  with  a  slightly  higher  variation  for  plastic  strains. 
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Figure  1  GM  Biaxial  Strain-Rate  Machine 
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Figure  2  Schematic  Diagram  of  GM  Biaxial 
Strain-Rate  Machine 
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Axial  load  to  the  specimen  is  provided  by  a  12.5  inch  tore 
cylinder  with  a  piston  constructed  of  titanium,  for  high 
strength- to-mass  ratio,  and  designed  to  operate  in  one  direc¬ 
tion  only.  The  piston  is  braced  between  the  stem  and  the 
flange  to  increase  its  stiffness,  the  stem  being  hollow  to 
allow  a  means  of  egress  for  instrumentation  cables.  Because 
of  the  unidirectional  movement  of  the  piston,  the  axial  com¬ 
pression  is  obtained  on  the  lower  end  of  the  piston  and  axial 
tension  on  the  upper  end.  Either  end  of  the  piston  rod  may 
be  attached  to  the  specimen  for  axial  loading,  with  a  variable 
hydraulic  damper  assembly  connected  on  the  other  end  to  damp 
undesirable  piston  r scillations . 

A  six  inch  bore  cylinder  is  used  to  load  an  auxiliary  piston 
which,  in  turn,  pressurizes  a  fluid  for  circumferential  load¬ 
ing  of  the  specimen.  The  specimen  is  pressurized  internally 
for  circumferential  tension,  and  externally  for  circumferen¬ 
tial  compression.  The  six  inch  pistons  are  not  braced  and  may 
be  operated  in  either  direction.  However,  because  of  the  uni¬ 
directional  nature  of  the  axial  loading  piston,  two  six  inch 
cylinder  assemblies  are  utilized,  one  for  circumferential 
loading  during  axial  tension  {upper  half)  and  one  for  circum¬ 
ferential  loading  during  axial  compression  (lower  half) . 

Fast  acting  valves  are  connected  to  the  front  reservoir  of 
each  cylinder  to  provide  a  means  of  evacuating  the  reservoir 
in  minimum  time.  Between  the  reservoir  and  the  fast  acting 
valve  is  an  interchangeable  orifice  plate  which  meters  the 
flow  of  gas  from  the  front  reservoir,  and  thus  acts  as  the 
primary  control  of  the  piston  velocity. 

The  machine  frame  is  designed  to  resist  separating  forces  of 
3 

400  X  10  lbs.  by  using  prestressed  column  tie-rod  construction. 
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The  construction  is  essentially  that  of  two  frames  with  a 

common  middle  plate.  As  only  one  end  of  the  machine  is  used 

at  a  time,  threading  the  tie  rods  through  the  middle  plate 

allows  each  end  of  the  tie  rod  to  be  stressed  separately. 

Calculations  show  deflections  of  the  order  of  0.002  inch 

3 

under  maximum  loading  (340  X  10  lbs.). 

Specimen  packages  were  designed  and  constructed  for  biaxially 
testing  tubular  specimens  in  each  of  the  four  stress  quadrants 
in  the  circumferential  stress — axial  stress  plane.  A  package 
serves  to  support  and  seal  the  tubular  specimen,  to  provide  a 
closed  pressure  chamber,  and  to  provide  support  for  pressure 
transducers  and  strain  gage  seals.  Each  package  consists  of 
a  piston  rod  adapter,  load  cell,  tubular  specimen  with  adapter, 
base,  ring  piston,  actuating  pins,  base  plug,  and  pressure 
transducer. 


Figures  3  and  4  show  cross  sections  through  two  packages ,  which 
differ  primarily  in  the  means  of  gripping  the  specimen  due  to 
the  nature  of  loading.  The  package  in  Figure  4  is  photographed 
in  Figure  5.  Both  of  these  packages  are  for  circumferential 
tension  and  thus  show  a  core  rod  through  the  center  which  serves 
to  align  the  assembly  and  reduce  the  volume  of  pressurizing 
fluid  required  to  stress  the  specimen.  Reduction  of  fluid 
volume  reduces  allowable  fluid  compressibility,  thereby  reduc¬ 
ing  oil  piston  stroke  requirements.  The  housing  to  contain  the 
fluid  in  the  circumferential  compression  package  serves  the 
same  purpose  as  the  core  rod.  Note  that  the  nose  of  the  ring 
piston  has  a  land  that  cushions  the  end  of  its  stroke  by  meter¬ 
ing  fluid  flow  between  the  base  and  the  core  rod.  The  metering 
arrangement  reduces  the  possibility  of  catastrophic  piston  im¬ 
pact  if  premature  failure  of  the  specimen  occurs  during  a  test. 
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Figure  5  Photograph  of  a  Test  Package  for 
Metal  Specimens 
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The  machine  is  operated  by  charging  each  cylinder  with  gas 
(air,  nitrogen,  or  helium)  introduced  at  equal  pressure  into 
a  large  reservoir  in  back  of  the  piston  and  into  a  small  re¬ 
servoir  in  front  of  the  piston.  The  piston  moves  forward 
when  the  small  reservoir  is  exhausted  by  flow  through  an 
orifice.  Exhaust  is  initiated  through  the  opening  of  the 
downstream  fast  acting  valve.  Piston  velocity,  and  hence  the 
rate  of  loading,  is  controlled  by  the  type  of  pressure  of  the 
working  gas,  the  orifice  size,  and  to  some  extent,  the  speci¬ 
men. 

A  computer  program  of  the  piston  motion,  described  in  detail 
in  Reference  1,  is  used  to  determine  the  initial  conditions 
of  the  driving  gas,  the  gas  pressure  required,  orifice  size, 
and  amount  of  viscous  damping  necessary  to  biaxially  load  a 
specimen  at  a  predetermined  stress  ratio.  Briefly,  the  pro¬ 
cedure  is  as  follows.  Uniaxial  stress  data  for  the  specimen 
material  is  used  to  produce  an  effective  stress-strain  curve 
(from  the  octahedral  stress-strain  relations)  for  the  material, 
with  the  uniaxial  stress  tensile  and  compressive  fracture 
stresses  as  boundary  conditions,  if  appropriate.  Through  a 
short  iteration  program  the  individual  stress-strain  data  for 
a  given  stress  ratio  are  obtained  for  each  direction  of  load¬ 
ing.  The  equations  of  continuous  curves  fitted  to  these  data 
are  used  i.i  the  specimen  reaction  calculations  in  the  piston 
motion  study.  ^  In  this  program,  the  input  parameters  of 
gas  pressure,  orifice  size,  and  damping  coefficients  are 
iterated  upon  until  the  desired  constant  strain  rate  is  ob¬ 
tained.  These  final  parameters  are  then  the  machine  settings 
used  in  the  actual  experiment.  Through  this  programming, 
efficient  use  is  made  of  machine  time  and  specimens  by  vir¬ 
tually  eliminating  trial  and  error  adjustment  in  the  system. 
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Upon  firing,  the  test  specimen  tube  is  subjected  to  axial  load 
by  displacement  of  the  main  piston  rod  of  the  Biaxial  Machine. 
Displacement  of  the  six  inch  cylinder  piston  rods  exerts  a 
force  on  the  actuating  pins  of  the  package,  and  these  in  turn 
drive  the  ring  piston  in  the  package  which  compresses  the 
fluid  in  the  oil  chamber.  The  pressurized  fluid  develops  the 
circumferential  stresses  in  the  tubular  walls. 

Instrumentation  for  the  Biaxial  Machine  consists  of  an  axial 
load  cell  and  a  pressure  transducer  for  loading  measurement, 
and  axially  and  circumferentially  oriented  strain  gages  to 
record  the  strain  on  both  the  inside  and  outside  surfaces  of 
the  specimen.  The  output  from  each  transducer  is  continuously 
monitored  in  analog  form  on  a  fourteen  channel,  wide  band 
(400  KHz),  frequency  modulated  tape  recorder.  Between  the 
tape  recorder  and  the  transducer  is  a  system  of  bridge  and 
balance  units,  signal  amplifiers  and  shunt  calibration  resis¬ 
tors  . 

An  IBM  1800  computer  is  used  as  a  central  processor  to  reduce 
the  measured  data  after  it  has  been  converted  from  analog  to 
digital  form.  This  conversion  technique,  coupled  directly  to 
the  computer,  omits  all  intermediate  steps  required  for  manual 
digitizing,  reducing  the  cumulative  system  error  to  approxi¬ 
mately  ±  0.5%  of  full  scale.  The  particular  processor  used 

3 

has  a  16  X  10  ,  sixteen  bit  word  storage  as  well  as  an  associ- 

g 

ated  10  word  two  disk  storage  unit.  This  flexible  system, 
shown  schematically  in  Figure  6,  has  the  capability  to  digit¬ 
ize,  store,  manipulate  with  the  proper  scale  factors,  print 
and  plot  the  various  reduced  data.  More  complete  details  are 
presented  in  References  1  and  2. 
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Figure  6  Block  Diagram  of  Data  Measurement  and 
Reduction  for  Biaxial  Machine 

SPECIMENS 

All  specimens  reported  were  machined  from  three  inch  rod  stock 
2024-T3  aluminum  obtained  from  the  same  mill  run.  A  complete 
description  of  the  material  appears  in  Appendix  A.  Separate 
specimen  configurations  were  used  for  obtaining  data  in  each  of 
the  four  stress  quadrants.  The  same  configurations,  however, 
were  used  for  both  the  virgin  and  prestrained  material  in  each 
quadrant,  with  the  exception  of  the  first,  or  tension-tension, 
quadrant  where  shorter  specimens  were  used  for  the  prestrained 
material.  In  machining  all  specimens,  diameters  were  held  to 
within  ±  0.0005  inch,  concentricity  of  inside  outside  diameter 
was  held  to  0.001  T.I.R.,  and  the  lengths  were  held  to  within 
0.03  inch  with  parallel  ends. 

For  the  virgin  material  in  the  tension-tension  quadrant,  the 
specimens  were  10.5  inches  in  length,  with  a  1.5  inch  I.D.  and 
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a  wall  thickness  of  0.050  inches.  The  ends  were  threaded  ex¬ 
ternally  to  provide  a  grip  for  transmitting  the  axial  load. 

This  left  an  effective  gauge  length  of  approximately  seven  inches 
To  avoid  column  buckling  in  the  compressive  prestrain,  an  explan¬ 
ation  of  which  follows,  billets  9.250  inches  in  length  had  to  be 
used,  necessitating  a  shorter  specimen  of  8.0  inches  in  length 
with  a  5.0  inch  effective  gauge  length.  All  specimens  for  the 
second,  or  axial  compression-circumferential  tension  quadrant, 
were  2,5  inch  O.D.,  2.25  inch  I.D.  tubes  5.0  inches  in  length. 

No  threads  were  necessary  due  to  the  compressive  axial  load. 

The  same  specimen  dimensions  are  used  for  the  compression-com¬ 
pression  or  third  quadrant.  The  size  of  the  fourth  quadrant 
specimen  was  the  same;  however,  this  specimen  was  threaded  in¬ 
ternally  to  provide  an  axial  grip,  and  then  had  the  outer  dia¬ 
meter  reduced  to  2.4  inches  over  a  central  3  inch  section  of  the 
cylinder  to  reduce  the  strength  requirements  on  the  threaded  ends 
Figure  7  shows  two  of  the  specimens,  both  tension-tension 
and  compression- tension,  to  add  visual  representation  to  the 
above  discussion. 

Additional  tests  were  made  using  an  Instron  uniaxial  testing 
machine  to  define  the  yield  surface  intercepts  on  the  axes  of 
the  stress  plane.  For  the  axial  stress  axis  the  specimens 
were  cut  parallel  to  the  axis  of  the  aluminum  rod,  while  for 
the  radial  stress  axis  the  specimens  were  cut  in  the  trans¬ 
verse  direction,  perpendicular  to  the  rod's  longitudinal  axis. 
Compression  tests  were  made  using  a  specimen  0.375  inches  in 
diameter  and  0.50  inches  in  length.  Threaded  specimens  with 
an  overall  length  of  2.625  inches,  gauge  diameter  of  0.25 
inches,  and  a  gauge  length  of  1.0  inches  were  used  for  tensile 
tests . 
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Figure  7  Second,  Third,  and  Fourth  (see  text) 

Stress  Quadrant  Specimens  on  the  Left, 

First  Quadrant  Specimen  to  the  Right 

All  prestraining  of  the  specimens  was  done  by  upsetting  billets 
cut  from  the  3.0  inch  diameter  rod  in  a  2000  ton  mechanical 
press  equipped  with  parallel  flat  platens.  "Barreling"  of  the 
specimens  was  minimized  by  coating  the  platens  with  a  Teflon 
based  metalworking  lubricant.  Variation  of  diameter  along  the 
billet  axis  did  not  exceed  0.015  inch  with  the  maximum  variation 
occurring  near  the  ends,  which  were  subsequently  removed  during 
machining.  For  the  second,  third,  and  fourth  quadrant  prestrained 
specimens,  6.0  inch  long  billets  were  given  an  axial  compressive 
prestrain  of  8.2%  by  subjecting  the  ends  to  an  axial  load  of 
225  tons.  Billets  9.25  inches  in  length  were  compressed  to  8.5 
inches  to  give  an  axial  compressive  '.strain  of  8.1%  for  the  pre¬ 
strain  specimens  used  in  the  tension- bension  quadrant.  Following 
the  upsetting  process,  specimens  were  machined  according  to  the 
configurations  previously  outlined. 
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DATA  REDUCTION  AND  RESULTS 


The  tubular  specimens  with  d/t  <  20  were  strain  gauged  on  the 
inner  and  outer  surfaces  with  90°  rosettes  to  measure  both  axial 
and  circumferential  strains.  Tubes  having  d/t  s  20,  those  used 
in  the  tension- tens ion  quadrant,  were  gauged  on  the  outside  sur¬ 
face  only,  because  variation  of  strains  through  the  thickness 
is  small.  Data  was  recorded  from  the  strain  gauges,  the  cali¬ 
brated  axial  load  cell,  and  the  fluid  pressure  monitoring  trans¬ 
ducer.  These  data  were  later  played  back  through  the  analog  to 
digital  converter  into  the  computer. 

A  complete  description  of  the  data  analysis  technique  is  given 
in  Reference  3,  however,  a  short  synopsis  follows  here  The 
data  interpretation  is  based  on  the  assumptions  of  homogeneous 
materials  following  conventional  elasticity  and  plasticity. 

The  subscripts  a,  t,  and  r  refer  to  axial,  circumferential, 
and  radial  respectively;  all  stresses  are  true  stress  values, 
and  all  strains  are  natural  or  logarithmic  strains  which  are 
given  by : 

e  =  ln{e  +  1)  e  =  engineering  strain 

e  =  natural  strain  (I-l) 

For  d/t  £  20,  thin  walled  tube  equilibrium  relations  were  used 

to  obtain  a  ,  with  a  being  calculated  as  F/A  (force/ins tantan- 

eous  area).  In  the  case  of  d/t  <  20,  variation  in  the  stresses 

through  the  thickness  need  be  considered  and  Lame's  elasticity 

relations  were  used  to  calculate  a  and  o J_.  A  discussion  of 

r  t 

the  error  involved  in  calculating  stresses  is  given  in  Refer¬ 
ence  3.  The  stresses  were  used  to  calculate  the  octahedrax 
shear  stress,  which  is: 
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To  *  1/3  *  (aa'0t)2  +  (W2  +  (or-aa)2 


(1-2) 


For  the  thicker  walled  cylinders,  d/t  <  20,  the  radial  strain 
at  a  given  radius,  r,  was  computed  as: 


2(e.  -  e.  )  r  2r.2  .  , 

t-  t.  o  l  i  =  inner  radius  values 

e  =  £•  +  - /-■  5-  —  -  a -  „  ..  ,  d-3) 

t  -  2,  2  _  2.  °=  outer  radius  values 

1  £i  o 


which  is  strictly  an  elastic  calculation,  although  it  is  a  good 

(3) 

approximation  (errors  of  about  10%)  when  extended  into  the 
plastic  region  for  a  few  percent  plastic  strain.  On  the  thin 
walled  specimens,  where  e.  and  e  were  not  measured  on  the  in- 

u  ci 

side  surface,  er  was  calculated  from: 


l-2v 

E 


(a  +  a  +  a  )  -  (e.  +  e  ) 
a  t  r  t  a 


(1-4) 


where  v  =  Poisson's  Ratio 
E  =  Young’s  Modulus 


This  relationship  assumes  that: 

(A)  The  strains  are  separable  into  elastic 
and  plastic  parts. 

(B)  The  sum  of  the  plastic  strains  is  zero. 

(C)  Elastic  strains  are  given  by  Hooke's  law. 

(D)  Strains  and  displacements  are  small. 
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The  strains  were  used  to  calculate  octahedral  shear  strain, 
which  in  terms  of  the  principal  normal  strain  is: 

Yq  =  2/3  V  <ea"et)2  +  (£t”er)2  +  W  (I-5) 

The  pressures,  loads,  stresses  and  strains  were  plotted  out 
automatically  after  computation,  with  the  variables  as  the 
ordinate  and  time  as  the  abscissa.  In  addition,  plots  were 
made  of  one  variable  against  another.  The  yield  surface 
was  generated  by  the  use  of  an  offset  technique  on  the  oct¬ 
ahedral  shear  stress-shear  strain  plot.  Yield  was  defined 
as  the  intercept  on  the  stress-strain  curve  of  a  straight 
line,  offset  on  the  strain  axis  at  0.25%  strain  (equivalent 
to  a  0.2%  strain  on  a  uniaxial  test),  and  drawn  parallel  to 
the  linear  portion  of  the  cu~ve.  On  the  octahedral  curves, 
the  slope  of  this  elastic  portion  is  equal  to  the  shear  modu¬ 
lus,  G,  defined  as: 

G  =  E/ [2(1  +  v)]  (1-6) 

tq  =  Gyq  in  the  elastic  range.  (1-7) 

The  shear  stress  value  determined  in  this  manner  was  referred 
to  the  shear  stress  vs.  time  plot  to  find  the  time  of  yielding. 
This  time  was  then  used  to  find  the  value  of  all  other  vari¬ 
ables  at  yield. 

Plots  were  made  of  the  circumferential  stress  vs.  the  axial 
stress  to  determine  the  yield  locus.  A  least  squares,  ellipti¬ 
cal  fit  of  the  data  was  made  and  the  results  are  shown  in  Fig¬ 
ure  8  for  the  virgin  material  and  in  Figure  9  for  the  prestrained 
aluminum.  Note  that  the  virgin  data  ellipse  closely  approximates 
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would  equal  a  negative  one,  the  von  Mises  yield  criterion, 

T 

where  Al  and  A2  would  equal  one,  A^  would  be  a  negative  one, 

J  2 

A4  and  A5  would  be  zero  and  A6  would  equal  a  negative  Y  ,  where 
Y  is  yield.  The  two  constants  A4  and  A5  differ  most  from  the 
von  Mises  criterion,  which  has  the  effect  of  translating  the 
ellipse  center  slightly  away  from  the  origin.  For  the  pre¬ 
strained  data,  the  constants  differ  significantly  from  von  Mises, 
particularly  A3,  which  changes  the  elliptical  shape,  and  A4  and 
A5,  which  translate  the  ellipse.  This  is  to  be  expected  as  pre¬ 
straining  produces  anisotropy  in  aluminum.  All  data  points  re¬ 
present  quasistatic  behavior,  the  testing  having  been  carried 

-4 

out  at  strain  rates  approximately  10  /second.  Figure  8  and  9 
will  be  discussed  subsequently. 

Equation:  Al  0a2  +  A? 9 2  +  +  A4(Ja  +  A50{  +  A6  •  0 


Figure  8  Yield  Surface  for  2024-T3  Aluminum 
in  As-Received  Condition  (Virgin) 
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Equation:  A1  of  +  A20,2  +  A3  0,0,  +  A40,  +  A50j  +  A6  •  0 
0^  and  O',  are  in  KS I 

WhprA- 


Figure  9  Yield  Surface  for  2024 -T3  Aluminum 

Pre-Strained  by  8%  Axial  Compression 

To  incorporate  radial  stresses,  which  are  an  order  of  magnitude 
lower  than  the  others,  plots  are  made  on  the  it  plane,  or  the 
plane  showing  the  deviatoric  components  of  stress  only.  For 
isotropic  materials  following  the  von  Mises  yield  criterion, 
the  three  dimensional  yield  surface  is  a  cylinder,  and  all 
yield  points,  regardless  of  the  hydrostatic  component  of  stress, 
project  back  onto  the  tt  plane  in  a  circle,  its  center  being  at 
the  origin  of  the  three  Cartesian  stress  coordinates.  Since 
the  virgin  yield  data  does  not  vary  significantly  from  the 
von  Mises  criterion,  a  circle  has  been  fit  to  the  data,  as 
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shown  in  Figure  10.  In  the  case  of  an  anisotropic  material, 

the  yield  surface  is  no  longer  a  cylinder,  but  takes  on  a  more 

(4) 

random  shape,  as  shown  by  Hsu.  The  deviatoric  components  of 
stress  vary  with  the  hydrostatic  component,  therefore  the  yield 
locus  for  each  hydrostatic  pressure  projects  back  onto  a  unique 
curve  on  the  ir  plane.  Since  the  hydrostatic  pressures  for  each 
of  the  prestrained  specimen  yield  points  are  different,  projec¬ 
tion  back  to  the  it  plane  produces  a  series  of  random  points  not 
of  the  yield. 


Figure  10  7r-Plane  Yield  Surface  for  2024-T3 

Aluminum  in  As-Received  Condition  (Virgin) 


A  note  might  be  made  of  some  of  the  difficulties  experienced  in 
obtaining  the  data.  Buckling  became  a  severe  problem  in  the 
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third  or  compression-compression  quadrant,  with  the  specimens 
becoming  unstable  before  bulk  yielding  was  reached.  One  point 
was  obtained  for  the  virgin  material,  and  none  for  the  pre¬ 
strained.  Increasing  the  wall  thickness  to  mean  radius  ratio 
of  the  specimen  would  have  allowed  higher  stresses  before  buckl¬ 
ing  occurred,  however,  thicker  specimens  have  greater  variation 
of  stresses  and  strains  across  the  wall  and  reduce  the  resolu¬ 
tion  of  yield  and  thus  are  undesirable.  To  insure  a  valid  yield 
surface,  additional  points  were  obtained  in  the  first  and  second 
quadrants  and  an  existing  specimen  package  was  modified  to  ob¬ 
tain  data  in  the  fourth  or  axial  tension-tangential  compression 
quadrant. 


DISCUSSION 

Mair  and  Pugh^  present  four  possible  forms  of  subsequent  yield 
surfaces,  shown  in  Figure  11,  due  to  different  theories  of  work 
hardening.  One  of  these  is  the  isotropic  hardening  theory,  which 
assumes  that  the  subsequent  yield  surface  is  a  uniform  expansion 
of  the  initial  surface  about  the  same  center  point.  Note  that 
this  theory  allows  for  cross-effects,  or  cross  hardening,  wherein 
prestraining  in  one  direction  produces  strain  hardening  in  some 
other  direction.  This  is  exemplified  by  a  change  in  the  mag¬ 
nitude  of  a  vector  from  the  geometrical  center  of  the  yield  sur¬ 
face  to  the  surface  itself  following  prestrain  in  another, 
generally  orthogonal  direction.  Examples  of  cross-effects  are 
•  hown  in  Figure  12.  The  isotropic  strain  hardening  theory  does 
not  allow  for  the  Bauschinger  effect,  e.g.  the  reduction  in 
yield  stress  following  stress  reversal  (see  Figure  12) .  Another 
hardening  theory,  attributed  to  Batdorf  and  Budiansky,  is  the 
slip  or  independent  loading  surface  theory,  which  states  that  the 
subsequent  yield  surface  is  the  virgin  yield  surface  enlarged  to 
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include  the  prestrain  point  with  the  minimum  amount  of  area. 

It  leads  to  a  corner  on  the  yield  surface  and  does  not  allow 
for  cross-effects  or  fully  for  the  Bauschinger  effect.  A  third 
theory  is  the  kinematic  hardening  theory  which  assumes  that  the 
subsequent  yield  surface  is  formed  by  rigidly  translating  the 
virgin  surface,  without  expansion,  in  the  direction  of  the  pre¬ 
load.  There  is  allowance  for  the  Bauschinger  effect  but  not 
for  cross-effects.  The  final  theory  is  one  that  Mair  and  Pugh 
attribute  to  Hodge,  wherein  the  subsequent  yield  surface  is 
explained  by  a  combination  of  isotropic  and  kinematic  hardening. 
Both  cross-effects  and  the  Bauschinger  effect  are  covered  by 
this  general  theory,  which  allows  expansion  and  translation. 


STRESS  PLANE 


HARDENING  THEORY 
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LOADING  SURFACE 
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COMBINED  OR  GENERAL 
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Virgin  Yield  Surface _ 
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Figure  11 


Comparison  of  Hardening  Laws ,  or 
Subsequent  Yield  Theories 
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stress  plot.  They  saw  a  high  curvature  near  the  shear  stress 
axis  which  they  felt  was  evidence  of  a  corner.  Ivey^  in  an 
investigation  of  several  aluminum  alloys  obtained  results  simi¬ 
lar  to  Naghdi,  et  al.,  except  there  was  no  indication  of  a 
corner.  Both  investigations  found  the  yield  surfaces  to  be 
convex  although  not  necessarily  enclosing  the  stress  origin, 
and  both  found  the  initial  yield  surface  to  be  closely  approxi¬ 
mated  by  the  von  Mises  yield  theory.  From  these  it  would  ap¬ 
pear  that  kinematic  hardening  was  the  best  theory,  but  Szcze- 
pinski  and  Miastkowski v  tried  to  fit  the  subsequent  yield 
of  an  aluminum  alloy  with  a  pure  kinematic  hardening  model  and 
found  great  deviations  for  large  plastic  prestrains.  Working 
with  rectangular  specimens,  they  also  found  that  if  they  pre¬ 
strained  to  various  degrees  in  one  direction  and  then  loaded  to 
a  larger  strain  in  the  orthogonal  direction,  the  yield  surface 
was  largely  a  function  of  the  last  leg  of  the  prestrain  path, 
the  material  "forgetting"  its  previous  history.  Among  the 
studies  that  showed  stronger  evidence  of  corners  was  that  of 

(9) 

Bertsch  and  Findley  who  experimentally  were  able  to  produce 
"rounded  corners"  by  subsequent  yielding  under  approximately 
equal  biaxial  stress.  Shiratori  and  Ikegami^0^  working  with 
brass  found  that  corners  were  evident  for  proportional  stress 
loading  paths  but  were  not  observed  when  the  loading  path  was 
varied.  Both  investigations  found  that  corners  were  not  pre¬ 
sent  under  uniaxial  stress  prestrain  conditions. 

Aside  from  the  varied  materials  tested  by  the  investigators, 
the  major  cause  of  the  differences  in  the  reported  data  is  the 
assortment  of  yield  definitions  employed.  Mair  and  Pugh  indi¬ 
cate  this  in  discussing  why  tney  saw  cross-effects  in  copper 
while  Naghdi,  et  al.  and  Ivey  did  not.  The  latter  studies  de¬ 
fine  yield  as  the  first  departure  from  linearity,  requiring 
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extremely  sensitive  transducers,  while  Mair  and  Pugh  used  both 
a  "backward  extrapolation"  method  and  a  "proof  strain"  method. 

The  backward  extrapolation  method  extrapolates  the  stress-strain 
curve  back  to  the  intersection  of  the  elastic  and  plastic  tan¬ 
gent  lines,  while  a  proof  strain  is  a  pre-selected  offset  plastic 
strain.  Both  the  backward  extrapolation  method  and  the  proof 
strain  method  require  that  the  specimen  be  loaded  into  the  plas¬ 
tic  range,  producing  a  yield  stress  greater  than  the  linear 
limit  value.  The  definition  of  yield,  even  the  amount  of  off¬ 
set  in  a  proof  method,  plays  a  significant  part  in  whether  cor¬ 
ners  will  be  seen.  ^  In  the  investigations  above  where  proof 
strains  or  extrapolation  methods  were  used,  the  amount  of  plas¬ 
tic  strain  to  yield  was  approximately  .01-. 03%,  an  order  of  mag¬ 
nitude  lower  than  the  .2%  generally  employed  in  engineering  an¬ 
alysis  and  used  in  this  study  for  the  generation  of  design  data. 

Returning  to  Figures  8  and  10,  which  represent  the  virgin  yield 
surface  of  the  2024-T3  aluminum,  note  that  the  data  could  be 
well  represented  by  the  von  Mises '  yield  criterion.  The  stress 
plane.  Figure  8,  shows  the  ellipse  centered  slightly  away  from 
the  stress  origin  in  the  direction  of  the  positive  stress  axis , 
and  the  major  axis  of  the  ellipse  makes  an  angle  of  approximately 
45°  with  the  two  stress  axes;  thus  the  variation  from  von  Mises 
is  small.  In  manufacturing  the  aluminum  rod  stock,  the  rod  is 
stretched  slightly  along  the  axial  axis  accounting  for  the  slight 
positive  axial  offset  of  the  virgin  material  yield  locus  origin. 
The  it  plane  representation.  Figure  10,  shows  even  less  evidence 
of  the  origin  being  off  center,  and  therefore,  a  circle,  its 
center  at  the  origin,  has  been  fitted  to  the  data. 

Figure  9  provides  evidence  of  the  anisotropy  produced  by  pre¬ 
straining.  The  center  of  the  ellipse  is  significantly  offset 


I 
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in  the  direction  of  prestrain  and  the  major  axis  has  rotated 
away  from  the  45°  line.  The  ratio  of  the  minor  to  major  dia¬ 
meters  also  increases  indicating  a  general  distortion  of  the 
shape  of  the  von  Mises  ellipse.  Neither  an  isotropic  nor  kin¬ 
ematic  hardening  theory  would  describe  the  yield  surface  of  the 
prestrained  aluminum. 

Figure  13  shows  the  two  yield  surfaces  plotted  on  the  same  stress 
plane  for  comparative  purposes.  Both  yield  surfaces  are  convex 
and  enclose  the  stress  origin.  Recalling  the  previous  discussion 
on  yield  criterion,  note  that  the  surface  expands,  transates  in 
the  direction  of  the  compressive  prestrain-  and  also  rotates 
slightly.  Cross-effects  do  exist,  which  might  be  expected  from 
the  definition  of  yield  used  here,  for  Mair  and  Pugh  show  that 

2QZ4-B51  Aluminum 


CT(KSI) 


60  80 


Figure  13  Comparison  of  Virgin  and  Prestrained 
Yield  -Sux  faces 
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the  proof  strain  yield  definition  Jeads  to  cross-effects,  where 
the  departure  from  linearity  definition  does  not  produce  them. 
Likewise,  corners  detectable  using  departure  from  linearity 
would  not  generally  be  discernible  using  the  techniques  incor¬ 
porated  in  this  study.  Evidence  of  the  Bauschinger  effect  is 
seen  in  the  substantial  difference  in  the  compressive  and  axial 
yield  stress  following  compressive  axial  prestrain. 

The  subsequent  yield  surface  of  this  material  is  best  repre¬ 
sented  by  the  general  hardening  theory  attributed  to  Hodge. 

This  can  be  seen  by  a  comparison  of  the  lower  part  of  Fig¬ 
ure  11  with  Figure  13.  The  general  theory  is  a  combination 
of  the  isotropic  and  kinematic  (segmentwise  linear  theory) 
hardening  theories,  which  Hodge  called  the  piecewj.se  linear 
isotropic  hardening  theory.  As  stated  previously  it  is  the 
only  one  that  accounts  for  both  cross-effects  and  the  Baus¬ 
chinger  effect,  however  to  account  for  additional  distortion, 
such  as  the  rotation  seen  here,  Mair  and  Pugh  felt  the  gen¬ 
eral  theory  had  to  be  expanded.  Again,  the  general  theory 
does  provide  a  sufficiently  accurate  presentation  of  the  sub¬ 
sequent  yield  phenomena  to  be  of  use  for  the  2024-T3  aluminum. 


INTRODUCTION  AND  TECHNIQUE  DESCRIPTION 

In  many  applications,  structural  materials  are  heated  rapidly 
and  either  simultaneously  or  after  some  short  period  are 
loaded  to  high  stress  levels.  In  this  study,  yield  and  flow 
stress  were  measured  for  two  aluminum  alloys,  6061-T6  and 
2014-T6  at  elevated  temperature  (see  Appendix  A  for  complete 
material  descriptions) .  The  efteuL  of  heating  rate  was  studied 
at  temperatures  from  400  to  700°F  where  uniaxial  stress  tests 
were  conducted  at  a  nearly  constant  strain  rate.  Prior  work^11^ 
led  to  expectations  of  increas  .'1  tensile  yield  strength  with 
decreased  heating  time  (or  increased  heating  rate) . 

The  tests  were  conducted  utilizing  the  General  Motors  Ultra- 

(12  13) 

High  Temperature  Medium  Strain-Rate  Machine  '  where  direct 
electrical  resistance  is  employed  to  heat  the  test  specimen  (Fig 
ure  14).  Nearly  linear  heating  rates  in  the  range  of  0.25  to 
2500°F/sec  can  be  obtained  by  controlling  the  current  through 
the  specimen.  Immediate:  _  (within  several  milliseconds)  after 
reaching  test  temperature,  the  specimens  were  loaded  in  uni¬ 
axial  stress  at  a  strain  rate  of  about  10/second.  This  strain 
rate  was  used  for  all  tests  in  order  to  keep  the  loading  time 
(0.005-0.010  sec)  significantly  smaller  than  the  shortest  heat¬ 
ing  time  (0.100-0.200  sec). 
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At  heating  rates  below  0.25°F/sec/  it  was  necessary  to  use  a 
variation  of  the  linear  heating  method  described  above.  These 
low  rates  required  heating  the  specimen  over  a  period  of  hours 
or  days.  For  these  rates,  the  specimens  were  first  heated  in 
a  radiant  oven  to  the  test  temperature,  maintained  at  this  tem¬ 
perature  for  the  desired  time,  removed  from  the  oven,  placed  in 
the  testing  machine,  reheated  at  25°F/sec,  and  tested  as  des¬ 
cribed  above. 


Positive 


Figure  14  Schematic  of  Ultra-High  Temperature 
Strain-Rate  Machine 
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Manual  control  of  the  heating  rate  and  machine  firing  proved 
satisfactory  at  rates  up  to  250°F/second.  At  higher  rates,  a 
simple  open  loop,  timed-logic  control  was  used.  Consistant 
heating  rates  and  test  temperature  was  achieved  after  the  pro¬ 
per  equipment  control  settings  were  empirically  determined. 

A 

Specimen  load  was  measured  from  the  output  of  semiconductor 
strain  gages  mounted  on  the  test  machine  load  cell,  which  were 
energized  and  monitored  by  an  AC  carrier  oscilloscope  plug-in 
module.  Specimen  strain  in  the  reduced  uniform  gauge  section 
was  determined  using  two  electro-optical  trackers  operated 
differentially.  Dynamic  load  and  gauge-section  elongation 
data  were  observed  on  a  dual  beam  oscilloscope  and  recorded 
photographically . 

Specimen  temperature  at  heating  rates  of  250°F/sec  and  below 
was  measured  with  a  0.003  inch  diameter  chrome 1-alumel  thermo¬ 
couple.  The  thermocouple  was  mounted  with  a  slight  spring 
pressure  contact  on  the  specimen  surface  at  the  longitudinal 
center  of  the  gauge  length.  The  spring  pressure  was  provided 
by  the  thermocouple  leads .  This  technique  along  with  the  small 
thermocouple  diameter  provided  adequate  response  at  heating 
rates  to  250°F/sec  and  minimized  thermal  gradients  in  the  speci¬ 
men  due  to  sensor  heating  sinking.  No  embedding  holes  or  bond¬ 
ing  agents,  which  would  create  local  stresses  in  the  specimen, 
were  required. 

Thermocouple  output  voltage  was  observed  on  a  chart  recorder 
and  tracked  against  a  pre-determined  graph  to  achieve  a  uni¬ 
form  heating  rate.  At  rates  above  250°F/sec  the  response 
limitations  of  the  thermocouple  and  chart  recorder  required  a 
change  in  technique  for  temperature  measurement.  A  four  mic¬ 
rosecond  time  response  infrared  radiometric  microscope,  stati- 
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cally  calibrated  with  a  thermocouple  reference  on  a  typical 
test  specimen,  was  used  and  observed  on  an  oscilloscope. 

The  extremely  low  electrical  resistance  and  high  thermal  con¬ 
ductivity  of  the  aluminum  alloy  required  minimizing  the  cross- 
sectional  area  of  the  test  specimen  to  obtain  rapid  and  uni¬ 
form  heating.  A  diameter  of  0.100  inch  was  selected  as  a 
practical  minimum  to  eliminate  possible  bending  of  the  speci¬ 
men  during  fabrication,  handling,  and  insertion  into  the  test¬ 
ing  machine.  In  addition,  this  diameter  was  significantly 
greater  than  the  mean  grain  size  of  the  alloys.  Machining  of 
the  specimens  of  the  standard  configuration  (Figure  15a)  with 
a  gauge-section  diameter  of  0.100  inch  was  impractical.  There¬ 
fore,  a  smaller  specimen  (Figure  15b)  was  designed  along  with 
re-usable  adapters  (Figure  15d)  to  fit  the  specimen  to  the  test 
machine  grips .  Specimen  blanks  were  cut  parallel  to  the  roll¬ 
ing  direction,  or  x-axis  (see  Appendix  A  -  Figure  A3) ,  of  the 
0.5  inch  2014-T6  plate  by  sawing  and  then  lathe- turned  to  a 
cylindrical  shape.  Specimens  (Figure  15b)  were  machined  from 
these  cylindrical  shapes,  and  from  the  6061-T6  rod  stock,  by 
precision  centerless  grinding.  Diameter  variations  along  the 
gauge  length  of  any  specimen  were  not  measurable  with  a  0.0001 
inch  resolution  micrometer,  and  varition  from  specimen  to  speci¬ 
men  was  consistantly  less  than  0.0005  inch.  Concentricity 
along  the  specimen  was  within  0.001  inch  TIR  from  end  to  end 
and  surface  finish  was  typically  two  to  three  microinches. 


Note:  Specimens  are 
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TEST  RESUT  V: 

Raw  test  data  in  the  form  of  oscilloscope  photographs  and 
micrometer  measurements  of  the  specimen  was  reduced  to  tabu¬ 
lar  and  graphic  form  with  the  aid  of  a  digital  computer.  The 
photographs  of  load  and  displacement  along  with  individual 
transducer  calibration  records  were  digitized  by  means  of  an 
optical  tele-reader.  This  data,  combined  with  the  specimen 
dimensions,  was  used  to  obtain  a  graphic  plot  of  engineering 
stress-strain.  The  uniaxial  stress  yield  was  obtained  from 
these  plots ,  wnere  yield  was  defined  as  the  first  observed 
deviation  from  linearity  in  the  elastic  portion  of  the  stress- 
strain  curve.  A  minimum  of  three  data  points  (and  usually 
five  or  more)  was  obtained  and  used  to  calculate  a  mean  value 
and  standard  deviation  value  for  each  heating  rate/test  tem¬ 
perature  combination.  Standard  deviation  values  were  computed 
according  to  the  formula: 


where 

a  =  standard  deviation 
x 

n  =  number  of  stress  values  taken 

th 

a.  -  deviation  of  1  stress  value  from  mean 
1  stress  value. 

This  parameter  was  chosen  as  reprerenting  a  conservative  esti¬ 
mate  of  the  probable  error  band  for  the  measured  stress  values 

(14) 

m  view  of  the  relat:vely  small  number  of  tests. 
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A  more  rigorous  analysis  of  the  data  would  include  modifica¬ 
tion  of  the  "standard  normal"  distribution  to  encompass  the 
extremely  small  (<<30)  number  of  samples  taken.  The  "student 
t  distribution"  may  be  used  as  such  a  modifier. 


TABLE  I 

STUDENT  t  DISTRIBUTION 
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6 

1.440 

1.943 

2.447 

3.143 

3.707 

7 

1.415 

1.895 

2.365 

2.998 

3.499 

8 

1.397 

1.860 

2.306 

2.896 
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2.262 
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v-fl 

t2  “  2 

Student  t  Distribution  =  f(t)  c  (1  +  — ) 

f(t)  =  standard  normal  distribution 

c  -  constant 

v  =  degree  of  freedom  of  t 

=  n-1 

n  =  number  of  samples 


This  distribution  modification  is  used  to  calculate  a  confidence 
interval,  I,  according  to: 


I 


VF 


where 


(II-l) 


X  =  mean  value  (stress) 
a  =  standard  deviation 

X 

n  =  number  of  samples 

t „  =  parameter  given  by  "student  t  distribution", 
ct 

The  probability,  ,  that  x,  the  yield  or  ultimate  stress,  will 
be  in  the  interval  I  for  similar  tests  is  given  by: 


P  =  1  -  2a 

C) 


(II-2) 
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Equation  (ll-l)  is  of  the  form: 

I  =  X  ±  Cax 
t 

where  C  =  —  modifies  the  standard  deviation  interval  to  in- 

\nr 

elude  the  statistical  effect  of  small  sample  number. 

Table  I  lists  numerical  values  for  the  "standard  t  distribution" 

as  a  function  of  v  and  a;  Figure  16  shows  the  standard  deviation 

coefficient  C  as  a  function  of  the  number  of  data  samples  for 

confidence  levels  between  80%  and  99%  (0.8  <  P  <  0.99). 
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Figure  16  Standard  Deviation  Coefficient,  C,  as 
a  Function  of  Number  of  Samples  and 
Desired  Confidence  Interval 
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Thus,  for  -•  .  gn  purposes,  any  required  confidence  interval  for 
the  experii.  ^ntal  data  may  be  expressed  as  the  product  of  the 
standard  deviation  and  the  appropriate  coefficient  taken  from 
Figure  16.  It  is  obvious  from  this  figure  that  for  n  a  3  the 
standard  deviation  represents  a  conservative  estimate  of  ex¬ 
perimental  error  to  a  confidence  of  80%  at  worst. 

Typical  stress-strain  results  for  the  2014-T6  alloy  at  500°F 
are  shown  in  Figure  17.  Yield  and  flow  stress  as  well  as 
elongation-to-fracture  are  all  highly  dependent  on  the  heating 
rate.  Figures  18  and  19  are  plots  of  yield  and  ultimate  stress 
as  a  function  of  heating  rate  from  400°F  to  700°F.  The  smooth 
curves  were  drawn  in  general  to  represent  the  mean  values  of 
the  data  points.  The  error  bands  shown  represent  standard  de¬ 
viation  from  the  mean.  Similar  results  for  the  6C61-T6  alloy 
are  shown  in  Figures  20  and  21,  where  yield  and  ultimate  stress 
are  plotted  against  heating  rate  at  500°F  and  700*^  test  tempera¬ 
tures  . 

Elongation-to-fracture,  defined  as  the  total  elongation  of  the 
gauge  section  divided  by  the  initial  gauge  section  length,  for 
the  two  alloys  after  heating  to  500 °F  at  various  heating  rates 
is  shown  in  Figure  22.  In  general,  large  scatter  was  observed 
primarily  due  to  the  effect  of  current  flowing  through  the 
specimen  during  deformation  which  tended  to  heat  up  the  material 
in  the  necked  region.  This  current  flow  had  a  negligible  effect 
on  yield  (and  likewise  ultimate)  stress  since  little  area  re¬ 
duction  occurred  at  these  points  and  since  the  time  to  reach 
these  conditions  was  very  short.  A  test  was  conducted  to  mea¬ 
sure  the  temperature  rise  during  necking  of  the  specimen  center 
for  the  most  rapid  heating  test  (2500°F/second) .  Results  showed 
a  temperature  rise  in  excess  of  100°F  occurring  at  this  rate 
between  ultimate  stress  and  fracture.  Smaller  temperature  rises 
occurred  at  lower  heating  rates. 
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Figure  18  2014-T6  Aluminum  Alloy,  Yield  Stress  vs. 

Heating  Rates  at  Various  Temperatures 
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Figure  19  2014-T6  Aluminum  Alloy,  Ultimate  Stress  vs 

Heating  Rates  at  Various  Temperatures 
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Figure  21  6061-T6  Aluminum  Alloy,  Ultimate  Stress  vs. 

Heating  Rates  at  Various  Temperatures 
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Some  concern  arose  when  it  was  observed  that  the  data  indicated 
the  material  to  be  most  sensitive  to  heating  rates  between 
0.025°F/sec  and  0.25°F/sec  where  the  method  of  heating  was 
changed,  as  explained  previously.  To  investigate  the  effect 
of  these  technique  differences,  control  specimens  were  heated 
radiantly  with  a  nearly  linear  temperature  rise  of  0.025°F/sec 
and  then  tested  as  described.  The  strength  of  the  control 
specimens  (heated  linearly  as  a  ramp  function)  was  about  2  ksi 
(^10%)  greater  than  the  specimens  heated  at  constant  elevated 
temperatures  (approximately  a  step  function) .  This  small  dis¬ 
crepancy  may  be  interpreted  as  arising  from  the  increased  time 
above  the  artificial  aging  temperature  (,v350oF)  for  the  step 
function  heating.  This  time  above  the  critical  temperature 
would  be  in  error  at  most  by  a  factor  of  two  for  the  700°F 
tests,  which  is  insignificant  when  compared  to  the  seven  orders 
of  magnitude  of  heating  rate  studied.  Although  this  technique 
difference  is  evident,  the  resultant  error  is  small  and  does  not 
affect  the  trend  of  the  data. 


DISCUSSION 


The  principle  process  used  to  strengthen  T6  alloys  of  aluminum 
involves  artificial  aging  at  approximately  350 °F  following  solu¬ 
tion  hardening.  This  artificial  aging  accelerates  precipitate 
hardening.  The  smaller  percentage  of  alloying  elements  would 
tend  to  make  the  6061  precipitation  rate  lower  than  the  higher 
alloyed  2014.  Reactions  of  these  two  alloys  to  various 

heat  treatments  are  shown  in  Figure  23  which  indicates  only  a 
small  difference  in  rales  of  precipitation  of  6061  and  2014 
alloys.  The  time  required  for  a  maximum  hardening  is  slightly 
greater  for  6061  at  any  given  temperature  than  for  the  2014 
alloy,  but  the  difference  is  not  great. 
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Figure  24  compares  yield  stress  at  temperature  for  various  heat¬ 
ing  rates  to  published  data  for  6061-T6  (no  cold  work)  held  at 
elevated  temperature  for  times  between  one  half  hour  and  10,000 
hours. The  strength  at  500°F  and  700°F  for  heating  rates  of 
0.0025°F/sec  and  0.025°F/sec  is  in  general  agreement  with  the 
published  curves,  which  is  meaningful  in  that  the  time  to  reach 
test  temperature  is  of  the  order  of  50  hours  and  5  hours  respec¬ 
tively  at  tnese  rates.  At  heating  rates  above  0.25°F/sec,  yield 
stress  at  temperature  is  significantly  higher  than  the  yield 
stress  of  material  held  at  elevated  temperature  for  30  minutes, 
the  difference  being  about  50%  at  500°F.  As  is  seen,  the  trends 
of  the  alloy  strength  at  elevated  temperature  for  all  heating 
rates  are  similar  to  the  published  values  for  material  held  at 
constant  elevated  temperature  for  prolonged  periods  of  time. 
Figure  25  shows  similar  results  for  the  2014-T6  alloy  compared 
to  curves  representative  to  all  aluminum  alloys  in  the  T6  temper. 


TEMPERATURE  ("Fl 


Figure  24  Normalized  Yield  Stress  vs.  Test  Temperature 
for  6061-T6  Aluminum  After  Various  Times  at 
Temperature 


47 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-70-12 


o 


Figure  25  Normalized  Yield  Stress  vs.  Test  Temperature 
for  2014-T6  Aluminum  After  Various  Times  at 
Temperature 

Figures  26  and  27  illustrate  the  effect  of  temperature  and  heat¬ 
ing  race  on  precipitation  observed  during  this  study.  A  speci- 

* 

men  from  each  heating  rate  group  at  test  temperatures  of  500 °F 
and  700CF  was  sectioned  parallel  to  the  cylindrical  axis,  me¬ 
chanically  polished,  and  chemically  etched  using  Kellers  etch^^ 
for  60-90  seconds  to  accentuate  grain  boundaries  and  precipitates. 
The  prepared  surface  was  examined  optically  and  with  an  electron 
microscope.  Photographs  3  through  7  of  Figure  26  show  the 
6061-T6  alloy  after  heating  of  500°F  at  rates  from  0.0025  to 
2500 °F/ second.  Above  0.25°F/sec  (Photos  5,  6,  7)  the  precipi¬ 
tates  are  dispersed  so  finely  that  they  are  not  visible.  At 


These  specimens  were  heated  to  temperature  at  nearly  con¬ 
stant  heating  rate,  and  cooled  to  room  temperature  in  air. 
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heating  rates  of  0.0025°F/sec  (Photo  3)  and  0.025°F/sec 
(Photo  4),  relatively  large  precipitates  appear  distributed 
uniformly  throughout  the  grains  and  along  grain  boundaries. 
Photographs  1  and  2  show  the  same  trend  at  700°F  with  the  pre¬ 
cipitate  forming  rather  gross  inclusions  at  a  heating  rate  of 
0.025°F/sec  Photos  8  through  12  are  scanning  electron  microscope 
pictures  of  the  same  specimens  showing  the  material  texture  be¬ 
coming  very  coarse  and  porous  for  the  longer  heating  times;  the 
voids  are  believed  to  be  precipitate  locations  with  the  precipi¬ 
tate  removed  during  etching.  Again,  little  change  is  observed 
at  heating  rates  above  0 . 25°F/second.  No  evidence  of  grain 
growth  or  recrystallization  is  observed  except  perhaps  for  the 
700°F  tests  at  a  0.025°F/sec  (Photo  1)  heating  rate.  This  is 
in  good  agreement  with  published  observations  of  grain  growth 
in  aluminum  alloys  generally  above  750°F  but  not  below  650°F. 
Little  or  no  variations  in  precipitation  size  or  density  was 

seen  between  material  before  and  after  deformation  as  has  been 

(11) 

observed  previously. 

Microscope  examinations  of  the  2014-T6  alloy.  Figure  27,  shows 
the  same  trend  toward  precipitation  coalescence  with  longer 
heating  times  although  it  is  not  as  clear  as  the  6061  photos. 

The  2014,  due  to  the  high  percentage  (M.4%)  of  alloying  cop¬ 
per,  contains  relatively  large  particles  of  high  copper  pre¬ 
cipitate  in  the  T6  condition  as  evidenced  by  the  large  white 
spots  in  the  photographs,  point  1,  Photo  11,  Figure  27.  These 
precipitates  are  seen  to  grow  in  size  and  number  as  heating 
time  increases.  The  elemental  compositions  of  the  precipitates 
was  measured  in  a  relative  way  using  secondary  X-ray  emissions 
in  the  scanning  electron  microscope.  Point  1,  in  Photos  7 
through  11  represents  high  copper  precipitates  (^  5:1  ratio 
over  surrounding  area).  Point  2,  Photo  8,  indicated  high  man¬ 
ganese  content,  as  did  point  2,  Photo  9,  while  point  3,  Photo 
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8  showed  relatively  high  concentration  of  manganese,  iron, 
copper  and  silicon.  The  other  points  indicated  showed  some 
traces  of  the  alloying  elements  but  were  predominantly  alum¬ 
inum.  X-ray  analysis  of  the  6061-T6  alloy  was  thwarted  by 
the  small  percentage  of  alloying  elements  (1.2%  max)  and  by 
the  closeness  of  the  characteristic  X-ray  energies  of  the 
alloying  elements — the  X-ray  spectrometer  could  not  resolve 
the  aluminum,  silicon,  and  magnesium  lines  although  some  in¬ 
dication  of  silicon  precipitation  was  seen  at  various  pre¬ 
cipitate  spots  as  a  broadening  of  the  aluminum  Ka  X-ray  line 
on  the  high  energy  side.  The  silicon  Ka  line  (1.8  ev)  is  just 
slightly  higher  than  the  aluminum  Ka  line  (1.5  ev) . 

The  times  at  temperature  required  to  soften  the  material  are 
similar  to  the  aging  curves  of  Figure  23.  This  would  suggest 
that  in  heating  times  shorter  them  about  15  minutes ,  little 
variation  should  be  seen  in  the  mechanical  properties  as  a 
function  of  heating  rate  or  time  at  temperature.  Therefore, 
the  mechanical  properties  at  higher  heating  rates  might  be 
assumed  equal  to  those  between  0.25  and  2500°F/second.  This 
assumption  would  simplify  material  modeling  in  high  heating 
rate  environments,  but  its  validity  is  dependent  on  precipi¬ 
tation  hardening  being  a  predominant  strengthening  mechanism 
and  reaction  rates  being  too  slow  to  proceed  very  far  when 
the  material  is  heated  rapidly  and  for  short  periods  of  time. 
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Figure  26  Metallography  Study  of  6061-T6  Alumiirun 


Figure  27  Metallography  Study  of  2014-T6  Aluminum 
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SECTION  III 

SHOCK  PROPAGATION  AND  FRACTURE  CHARACTERISTICS 
OF  2014-T6  ALUMINUM 


INTRODUCTION 

The  prediction  of  the  response  of  structures  to  impulsive  load¬ 
ing  has  been  studied  extensively  during  the  last  decade.  Out 
of  these  studies  sophisticated  computer  codes  have  been  devel- 

(18  i 9  j 

oped  which  may  be  used  to  assess  structural  response. v  ' 

These  codes  are  designed  to  model  both  the  initial  phase  of 
loading,  when  stress  waves  are  propagated  through  the  struc¬ 
ture,  and  later  time  response,  when  large  permanent  deforma¬ 
tions  result.  However,  to  use  these  computer  codes  effectively, 
input  parameters  must  be  determined  for  the  structural  materials 
of  interest. 

The  purpose  of  the  present  study  is  twofold;  (1)  to  provide 
these  input  parameters  for  an  important  structural  aluminum 
alloy,  2014-T6,  and  (2)  to  provide  basic  understanding  of  the 
shock  propagation  and  fracture  characteristics  of  this  alloy. 

The  following  often  serve  as  input  data  for  modeling  of  shock 
propagation  and  dynamic  frac  are  (quasi-static  inputs  not 
listed):  elastic  constants,  hugoniot  eqi  )n  of  state,  pre¬ 

cursor  amplitude,  .nd  dynamic  spall  strength.  Values  for 
some  or  all  of  these  quantities  are  usually  sufficient  for 
making  preliminary  code  calculations  of  material  response. 
However,  if  predicting  ea _ ly  time  shock  propagation  is  cf 
interest,  time  dependent  models  usually  must  be  developed,  and 
to  assess  their  validity  certain  checks  must  be  applied.  Some 
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readily  available  checks  are  elastic  precursor  decay,  compres- 

* 

sive  wave  profiles,  spall  wave  correlation  profiles  ,  and 
release  paths. 

In  concluding  this  study,  an  effort  was  made  to  correlate  the 
results  to  available  data  on  another  aluminum  alloy,  60S1-T6 
aluminum.  It  is  hoped  that  by  such  a  comparison,  users  of 
these  reported  results  on  2014-T6  might  extend  their  under¬ 
standing  beyond  the  scope  of  this  study  to  the  wealth  of  avail¬ 
able  data  on  6061-T6  aluminum. 

EXPERIMENTAL  TECHNIQUES 

Experimental  measurements  made  are  separable  into  the  following 
categories:  elastic  constants,  hugoniot  equation  of  state, 
wave  profile  studies,  spall  recovery  tests,  and  spall  wave  pro¬ 
file  measurements.  Three  different  laboratory  instruments  were 
used  in  making  these  measurements :  ultrasonic  instrumentation 
for  elastic  constants  measurements,  a  63.5  mm  bore  compressed 
air  gun  for  spall  recovery  studies,  and  a  102  mm  bore  compressed 
air  gun  for  the  hugoniot,  wave  profile  and  spall  wave  profile 
measurements . 

1.  Elastic  Constants  Measurements 

The  three  basic  measurements  made  to  determine  elastic  con¬ 
stants  are  density  and  dilatational  and  shear  wave  velocity. 
Density  was  determined  by  measuring  the  mass  and  volume  of 
several  cylindricaily  shaped  specimens. 

* 

Spall  wave  correlation  profiles  refer  to  free  surface 
velocity-time  profiles  recordable  at  a  specimen  free 
rear  surface  beneath  which  dynamic  fracture  occurs. 
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Wave  velocities  were  determined  by  measuring  wave  transit 
times  through  specimens  of  known  thickness.  A  pulse  superposi¬ 
tion  technique  was  used  for  these  transit  time  measurements . ^ 20 ^ 

4 

This  technique,  which  gives  good  accuracy  1  part  in  10  ) , 
was  also  used  for  determining  changes  in  transit  time  as  a  func¬ 
tion  of  temperature  and  hydrostatic  piessure.  Measurements 
were  made  over  a  temperature  range  from  25  to  200 °C,  and  a  pres¬ 
sure  range  from  0  to  8  kilobars. 

A  block  schematic  of  pulse  superposition  instrumentation  is 
shown  in  Figure  28.  The  theory  of  operation  is  as  follows. 

A  transducer,  either  longitudinal  or  shear  mode,  is  bonded  to 
one  side  of  a  disk  shaped  specimen  with  parallel,  flat,  polished 
surfaces.  A  short  radio  frequency  burst,  matched  in  frequency 
to  the  resonant  frequency  of  the  transducer,  is  applied  to  the 
transducer.  The  transducer,  which  is  both  a  transmitter  and 
receiver,  transmits  an  ultrasonic  wave  into  the  specimen.  This 
wave  reflects  off  the  back  surface  of  the  specimen  and  is  re¬ 
ceived  as  an  echo  by  the  transducer.  If  a  second  RF  burst  is 
applied  at  the  time  this  first  echo  is  received  and  a  third  at 
the  time  the  second  echo  is  received,  etc,  the  echo  amplitudes 
will  add.  In  order  to  determine  when  addition  occurs,  it  is 
necessary  to  interrupt  the  RF  bursts  to  be  able  to  observe  the 
decaying  echos.  Then  the  repetition  rate  of  the  RF  bursts, 
which  determines  the  time  between  pulses,  is  adjusted  to  give 
maximum  echo  amplitudes,  i.e.,  exact  superposition.  The  wave 
transit  time  through  the  specimen  is  then  one-half  the  reciprocal 
of  the  oscillator  frequency  (which  controls  the  RF  repetition 
rate) .  For  the  experimental  setup  used,  the  variable  frequency 
oscillator  is  a  General  Radio  Model  116 3-A,  the  pulse  generator 
is  a  Hewlett-Packard  Model  214-A,  the  RF  pulse  oscillator  is  an 
Arenberg  Model  PG-650-C,  the  impedance  matching  transformer  is 
an  Arenberg  Model  WB-100,  the  wide  band  amplifier  is  an  Arenberg 
Model  WA-600-D  and  the  frequency  counter  is  a  TSI  Model  PA-620. 
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Figure  28  Pulse  Superposition  Technique  Schematic 


The  pulse  superposition  technique  was  also  used  to  determine 
tne  temperature  and  pressure  dependence  of  the  elastic  wave 
velocities,  however,  changes  in  specimen  length  with  tempera¬ 
ture  and  pressure  must  be  incorporated  into  results  of  transit 
time  measurements  in  order  to  compute  correct  wave  velocities. 

The  temperature  dependence  of  elastic  wave  velocities  is  deter¬ 
mined  by  immersing  the  specimen,  with  transducer  attached, 
into  a  constant  temperature  bath.  A  Blue  M  Electric  Co.  Model 
60  bath  with  HTF-100  UCON  fluid  was  used  for  this  purpose. 
Generally,  the  upper  temperature  limit  is  determined  by  the 
transducer/specimen  bond  material  and  is  about  300 °C  for  dila- 
tational  waves  and  100 °C  for  shear  waves. 

Measurement  of  the  pressure  dependence  of  wave  velocities  was 
made  using  a  static  high  pressure  apparatus  illustrated  in 
Figure  29.  The  hydraulic  system  consists  of  two  parts;  a  low 
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pressure  portion  which  transmits  pressure  to  a  high-pressure 
portion  which  transmits  pressure  to  the  vessel  containing  the 
specimen.  The  usable  portion  of  the  pressure  vessel  is  2.5  cm 
di^r  er  by  15  cm  long  and  the  working  fluid  is  Octoil-S.  A 
pre-  ,.re  range  from  0  to  8  kilobars  was  used  for  the  measure¬ 
ments  made  on  2014-T6  aluminum.  Pressure  was  measured  with  a 
manganin  pressure  cell  moni'ored  with  a  resistance  measuring 
bridge.  Since  pressure  dependence  is  determined  under  isother¬ 
mal  conditions,  temperature  of  the  working  fluid  was  monitored. 

Propagation  direction  was  normal  to  the  12.7  mm  plate  stock 
for  these  wave  velocity  measurements. 
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Figure  29  Static  High-Pressure  Apparatus  Schematic 


57 


MANUFACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 

MSL-70-12 

2.  Spallation  Threshold  Measurements 

For  the  spallation  studies  specimens  were  recovered  after  shock 
loading  and  examined  to  determine  the  extent  of  dynamic  frac¬ 
ture.  Spallation  recovery  tests  were  conducted  using  the 
63.5  mm  compressed  air  gun  shown  in  Figure  30.  These  shots  were 
instrumented  to  record  only  impact  velocity.  The  velocity  mea¬ 
surement  system,  ending  2  cm  from  the  muzzle  of  the  gun,  in¬ 
cludes  a  light  source,  a  collimating  lens,  and  five  sets  of 
equally  spaced  slits  in  the  launch  tube.  The  slits  create  five 
narrow,  parallel  light  beams  crossing  the  launch  tube  and  emerg¬ 
ing  from  holes  on  the  opposite  side,  where  they  are  focused  by 
another  lens  onto  the  sensitive  element  of  a  photomultiplier 
tube.  An  opaque  projectile  passing  through  the  launch  tube 
successively  interrupts  these  light  beams ,  so  that  when  the 
output  of  the  photomultiplier  is  recorded  with  an  oscilloscope, 
projectile  velocity  can  be  determined  to  within  ±  1%. 
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The  impactor  and  target  are  shown  schematically  in  Figure  31. 

For  the  spallation  threshold  studies,  the  impactor  and  target 
were  both  2014-T6  aluminum  with  the  ratio  of  target  thickness 
to  impactor  thickness  'v-  1.8.  Two  impactor  thicknesses  were 
used,  0.6  and  1.6  ran*.  The  extent  of  dynamic  fracture  was  deter¬ 
mined  by  sectioning  recovered  specimens  across  a  diameter, 
polishing,  etching,  and  examining  optically  at  a  magnification 
of  25  to  100X.  Specimens  were  then  graded  according  to  the 
degree  of  fractures  to  complete  material  separation.  The  in¬ 
cipient  spall  threshold  was  defined  as  the  impact  velocity  cor¬ 
responding  to  cracking  over  50%  of  the  specimen  on  a  projected 
length  basis. 


63.5mm 
LAUNCH  TUBE 


Figure  31  Room  Temperature  Spall  Test  Schematic 
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3.  Hugoniot  Measurements 

Four  shots  were  fired  to  determine  the  hugoniot  of  2014-T6  alum¬ 
inum  to  50  kilobars.  Two  of  these  shots  were  a  "direct  impact" 
technique  where  an  aluminum  specimen  was  impacted  directly  upon 
a  quartz  gauge.  The  other  two  shots  employed  a  "buffered  direct 
impact"  technique  where  a  tungsten  carbide  disk  was  placed  in 
front  of  the  quartz  gauge  in  the  target.  These  two  techniques 
are  illustrated  in  Figures  32  and  33.  The  buffered  technique 
was  used  for  the  higher  stress  shots  ('v  35  and  50  kilobars  be¬ 
tween  aluminum  and  tungsten  carbide) . 


impactor 

(Unknown) 


Current-Time 


Figure  32  Schematic  of  Direct  Impact  Technique 
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Figure  34  102  nun  Compressed  Gas  Gun  Schematic 


Projectile  velocity  was  measured  by  a  pin  shorting  technique. 

As  the  impactor  neared  the  target,  five  pins  of  accurately 
known  distance  separation  were  successively  shorted.  The  short¬ 
ing  time  of  each  pin  was  read  out  directly  on  a  digital  counter , 
from  which,  impact  velocity,  V^.,  was  computed.  Impact  velocity 
measured  this  way  is  accurate  to  about  0.5%. 

Tilt  was  not  recorded  for  hugoniot  shots  but  was  measured  on 
wave  profile  shots.  Tilt  was  measured  by  means  of  a  shorting 
pin  technique.  Four  pins,  all  in  the  plane  of  the  target,  are 
shorted  during  impact  by  the  face  of  the  impactor  which  is 
grounded  by  a  pin  protruding  ahead  of  the  tilt  pins.  From  tilt 
pin  shorting  times,  recorded  by  digital  counters,  tilt  angle  at 
impact  is  computed.  Tilt  angles  did  not  enter  into  data  analysis 
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but  were  used  only  to  determine  planarity  of  impact  for  the 
purpose  of  screening  out  bad  shots.  Tilt  angles  less  than  one 
milliradian  were,  in  general,  judged  acceptable  for  this  study. 
In  practice,  tilt  angles  were  typically  0.5  milliradians  or 
less. 

Quartz  gauges  used  for  stress  measurements  were  1/4"  X  1  1/4" 
diameter  disks  with  1/2"  active  electrode  diameter.  These 
gauges,  obtained  from  Valpey  Corporation  and  characterized 
by  Graham^2"^  ,  produce  a  current  proportional  to  stress  at 
the  front  surface,  for  a  recording  time  of  one  microsecond. 

The  active  electrode  area  of  each  gauge  was  connected  to  three 
oscilloscopes  using  low-loss  foam  dielectric  instrumentation 
cables.  The  oscilloscopes  were  Tektronix,  Type  547  and 
Type  454,  with  fast  writing  rate  P-11  phosphor.  Quartz  gauge 
current-time  data  were  recorded  on  Polaroid  film  along  with 
time  calibration  marks  and  current  calibration  lines. 

4.  Wave  Profile  Measurements 

Wave  profile  measurements  were  made  using  a  "transmitted  wave" 
technique.  A  fused  quartz  (tungsten  carbide  for  two  shots) 
impactor  was  propelled  against  an  aluminum  target.  After  im¬ 
pact,  stress  waves  propagated  through  to  the  target  back  sur¬ 
face,  where  they  were  recorded  by  measuring  either  stress 
with  a  quartz  gauge,  or  interface  velocity  between  the  target 
and  fused  quartz.  Interface  velocity  measurements  were  made 
with  the  laser  velocity  interferometer  for  which  fused  quartz 

was  used  as  a  "window"  for  the  laser  beam.  It  has  been  char- 

(22) 

acterized  for  this  purpose  by  L.  Barker. 
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The  instrumentation  system  for  the  laser  velocity  interferometer 
is  shown  in  Figure  35.  This  system  consists  of  a  Perkin-Elmer 
Model  5800  Laser,  a  series  of  mirrors  directing  the  laser  beam 
towards  the  target ,  a  focusing  lens  to  focus  the  beam  at  the 
specimen-fused  quartz  interface  in  the  target,  a  silvered  dot 
at  this  interface  from  which  light  is  reflected,  a  series  of 
mirrors  to  direct  the  reflected  beam  back  to  an  optical  bench, 
a  beam  splitter  which  splits  the  reflected  light  (half  of  which 
goes  through  a  delay  leg  but  then  is  combined  later  with  the  re¬ 
flected  beam) ,  and  an  ITT  Model  F4034  photomultiplier  tube. 


Figure  35  Laser  Velocity  Interferometer  Schematic 
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The  output  of  the  photomultiplier  tube  is  connected  to  Tektronix 
Type  519  and  454  oscilloscopes.  A  fringe  pattern  is  developed 
by  recombining  the  laser  beam  after  the  delay  leg.  This  fringe 
pattern  will  change  as  the  interface ,  from  which  the  beam  is 
reflected,  moves.  The  number  of  fringes  seen  by  the  PM  Tube  per 
unit  time  is  proportional  to  acceleration  of  this  interface.  Th 
proportionality  constant  is  controlled  by  the  optical  length  of 
the  delay  leg.  A  target  and  projectile,  instrumented  for  the 
laser  velocity  interferometer,  is  shown  in  Figure  36. 
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Figure  36  Projectile  and  Velocity  Inter 
ferometer  Target 
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Two  shots  were  fired  using  quartz  instrumentation  and  six  using 
the  laser  velocity  interferometer  for  wave  profile  measurement. 
Compressive  wave  behavior,  release  wave  behavior,  and  wave  at¬ 
tenuation  in  2014-T6  aluminum  were  studied  wiT.h  these  shots. 

On  all  but  two  shots  the  impactor  was  a  fused  quartz  disk 

1.5  mm  thick).  A  tungsten  carbide  impactor  1  mm  thick) 
was  used  for  shots  9780  and  9781.  The  stress  wave  input  into 
aluminum  is  shown  schematically  in  Figure  37  for  these  two  types 
of  impactors.  Fused  quartz  inputs  a  square  stress  pulse  to 
the  target,  whereas  tungsten  carbide  inputs  a  stepped  stress 
pulse  as  the  target  is  unloaded  in  discrete  stress  steps  due 
to  unloading  waves  reverberating  back  and  forth  in  the  impac¬ 
tor.  This  "reverberating  buffer"  technique  has  been  described 

t  (23) 

by  Lysne. 
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Figure  37  Stress  Inputs  for  Wave  Profile  Studies 
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5.  Spallation  Wave  Profile 

Spall  wave  profile  measurements  were  made  on  the  102  mm  gun 
using  the  laser  velocity  interferometer.  The  configuration 
for  these  shots  is  shown  in  Figure  38.  Aluminum  was  impacted 
onto  aluminum  with  readout  of  free  surface  velocity  at  the 
rear  surface  of  the  target.  Shots  were  fired  below  and  above 
the  incipient  spall  threshold  for  both  loading  times  studied 
previously  on  the  63.5  mm  gun.  Impact  velocity  and  tilt  were 
recorded.  Specimens  were  recovered  and  .ramined  metallograph- 
ically  to  determine  the  extent  of  fracture. 


~.6mm  ~1. 07  mm 

or~L6mm  or~3.Q5mm 


Figure  38  Schematic  of  Spall  Wave  Profile  Shots 
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RESULTS  AND  DISCUSSION 


1.  Elastic  Constants 


List  of  Symbols 


B 

B 

C 


T 


P 

T 

a 

3 

Y 

P 


Adiabatic  bulk  modulus 

Isothermal  bulk  modu?  us 

Longitudinal  wave  velocity 

Heat  capacity  at  constant  pressure 

Shear  wave  velocity 

Hydrostatic  pressure 

Temperature 

Coefficient  of  linear  expansion 
Coefficient  of  volume  expansion 
Gruneisen  ratio 
Density 


Results  of  ultrasonic  and  density  measurements  on  2014-T6 
aluminum,  12.7  mm  stock,  are  listed  in  Table  II.  Pressure 


3Ct 


3Cr 


"T  w  WC 

derivatives  (■%=— )„  and  (-r£-~)m  have  been  corrected  for  length 
3P  T  3P  T  3C  ac 

changes  under  hydrostatic  pressure,  (~=)p  and  (^~)p  have 

been  corrected  to  account  for  thermal  expansion. 
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TABLE  II 

RESULTS  OF  ULTRASONIC 
MEASUREMENTS  ON  2014-T6  ALUMINUM 


UNITS 


UNCORRECTED 


CORRECTED 


ran 

psec 


_ 1_ 

kbar 


ran 


psec  kbar 


ran 

psec 


_1_ 

’C 


ran 

psec 


I 

°C 


0.0188 


0.0118 


-10.6  X  10 


-4 


-7.8  X  10 


-4 


0.0158 


0.0103 


-9.2  X  10 


-4 


-7.1  X  10 


-4 


UNITS 


ran/psec 


mm/psec 


g/cm 


3 


WAVE  VELOCITY  AT 
25°C  AND  1  ATM 


6.353 


3.199 


2.81 


Derivatives  evaluated  at 
25°C  and  1  atm. 


* 


Corrected  for  length  changes 
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Elastic  moduli  and  their  derivatives  may  be  computed  from 
quantities  listed  in  Table  II.  This  has  been  done  for  the 
bulk  and  shear  moduli.  The  following  equations  were  used 
for  these  calculations : 


fiS  =  P<CL  - 


(III-l) 


rT  -  ■  3 

1  +  BYT 


(III-2) 


a»S  (  3CT  4  ^Co 

(3P~)T  =  2p  )  CL  (3P“}T  "  3  Cs  (8P~) 


T)  +  bt 


(III-3) 


4  3*“"g 

(3T_,P  =  2p  CL  (3T~)P  3  Cs  (3T"“)P 


-  8BL 


(III-4) 


(3£\  =  (Mi)  +  (3B®  TX 

V3P  's  V3P  jT  '3T  ;P  BS 


(III-5) 


/ 3BS -  .3BS,  f9BS, 

'  3T  S  Ty  '3P  jT  l3T  ;p 


(III-6) 


‘§V  ji-^rcS^p -*«#>,) 

+  (ByT^)2  |(#)t  -  1  -  ^2  <!£>} 

<|£>P  -  £  «S£>P  -  {"«&-*") 


(III-7) 


(III-8) 
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m  rn  rp 

(321)  =  (331)  +  (321)  iX 

V9P  'S  l3P  'T  '3T  'p  -S 

D 


(III-9) 


,33!)  .  3!  (32!)  +  ,33l, 

v  3T  '  S  Ty  '3P  'T  V3T  ;p 


(III-10) 


G  =  Cgp  =  GS  =  GT 


(III -11) 


(—)  =  r 

3T  P  G 


l«s 


3C_ 

-  6 


(III-12) 


,3Gx  _  r  (  _2  .  1_ 

^  3P'  T  I  Cg  V3P  '  T  fiT 


(III-13) 


f3Gj 

1 3P'  S 


(3G)  +  ll  (36) 

'  3P;  T  DS  '3T*  P 

15 


(III-14) 


(36) 

V3TJS 


(32,  +  £  (32) 

' 3T' P  Ty  '  3P '  T 


(III-15) 


Values  for  the  above  quantities  are  given  in  Table  III,  all 
derivatives  listed  were  evaluated  at  T  =  20 °C  and  P  =  0  kilo- 
bars.  In  calculating  these  quantities,  values  given  below 
were  used  for  8  and  y.  The  volume  coefficient  of  expansion, 


8,  was  taken  as 


(24) 


8  =  3a  =  (64.2  +  0.039T)  X  iO_6/°C 


(III-16) 


The  Gruneisen  ratio,  y,  was  calculated  as 


y  =  1.81 


(III-17) 
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from  the  relation 


Y  = 

T  PC. 


(III-18) 


with  BS  taken  from  Table  III  and  CD  =  0.23  cal/g-°C. 


TABLE  III 


ELASTIC  CONSTANTS  FOR  2014-T6  AL 
(P  =  0,  T  =  25°C) 


PARAMETER 


751  kbar 


-0.207  kbar/°C 


5.68  kbar/'C 


-0.296  kbar/'C 


5.82  kbar/°C 


-0.147  kbar/'C 


2.99  kbar/'C 
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Several  different  equations  of  state  may  be  computed  from  the 
moduli  and  derivatives  of  moduli  listed  in  Table  III.  In 
the  discussion  that  follows,  isotherm,  adiabat  and  hugoniot 
calculations  are  made  from  the  Murnaghan  or  logarithmic  equa¬ 
tion  of  state,  for  which  pressure  along  an  isotherm  is  given 
by 


B 

P  =  T 
T  ,9BX, 

'  9P  '  T 


T 

V  ,9b\ 

(^)  '  9P  T  -  1 


(III-19) 


pressure  along  an  adiabat  is  given  by 

„S 


P  =  -51 

S  <*S°> 

'  9P 


(  V  ,  9B^.  ) 

9P  S  -  1  | 


(III-20) 


and  hydrostatic  pressure  along  the  hugoniot,  i.e.  shock  hydro- 
stat,  is  given  by 


p  =  P  +  IE 
H  S  12 


I 

+ 1 


a  -  f-)3 

o 


(III-21) 


Equations  of  state,  calculated  from  Eqs.  (III-19) ,  (III-20) 
and  (III-21)  are  presented  in  Figure  39  for  2014-T6  aluminum. 


For  pressures  up  to  50  kilobars  the  hugoniot  calculated  from 
Eq.  (III-21)  lies  within  0.2%  of  the  adiabat,  whereas  the  separa¬ 
tion  between  the  adiabat  and  isotherm  is  %  3%.  The  conversion 
to  particle  velocity  from  volume  was  made  using 
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2.  Spallation  Threshold 

Spall  fracture  characteristics  of  2014-T6  aluminum  are  sum¬ 
marized  in  Table  IV  and  shown  in  Figures  40,  41  and  42. 

Two  loading  times  were  studied  corresponding  to  impactor  thick¬ 
nesses  of  ^  0.6  mm  {0.19  ysec  pulse  duration)  and  ^  1.6  mm 
(0.5  ysec  pulse  duration).  Spall  fracture  characteristics  are 
shown  in  Figure  40  for  the  0.6  mm  impactor  series  and  in  Fig¬ 
ure  41  for  the  1.6  mm  impactor  series.  Spall  fractures  for 
each  follow  the  direction  of  grain  elongation  although  there 
is  a  tendency  for  cracks  to  link  up  at  right  angles  to  this 
direction.  Incipient  spall  thresholds  were  estimated  as 
0.216  mm/ysec  and  0.180  mm/ysec  respectively  for  these  two 
impactor  series. 

TABLE  IV 

SPALL  RESULTS  FOR  2014-T6  ALUMINUM 


Shot  No. 

Velocity 

(ram/usec) 

! 

Spall 

Projectile 

Thickness 

(mm) 

Target 

Thickness 

(mm) 

115 

.237 

o 

1.59 

3.04 

102 

.200 

© 

1.57 

3.06 

107 

.183 

0 

1.58 

3.03 

106 

.180 

® 

1.57 

3.04 

104 

.170 

e 

1.57 

3.05 

105 

.164 

0 

1.59 

3.05 

114 

.270 

• 

0.605 

1.067 

109 

.247 

© 

0.597 

1.072 

110 

.242 

9 

0.597 

1.069 

111 

.216 

® 

0.607 

1.069 

112 

.190 

O 

0.605 

1.072 

113 

.170 

o 

0.605 

1.067 

•  Complete  Separation 
©  Above  Incipient 
®  Incipient  Spall  Threshold 
Q  Below  Incipient 
O  No  Spall  Fracture  Visible  at  40X 


mm 
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Photomicrographs  at  40x 
From  Recovered  Specimens 
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113 


*  Spail  Threshold  Estimated  at  0. 21610. 005mm/*»sec 


impact  Velocity 


Q.  247mm/*/ sec 


0.242mm/*/ sec 


0. 216mm/*/ sec  * 


0.19ftnm/*zsec 


0. 170mm/*isec 


Figure  40  Metallographs  of  Spall  Fracture  in 
2014-T6  Aluminum,  0.6  mm  Impactor 
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Photomicrographs  at  40X 
From  Recovered  Specimens 
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Figure  41  Metallography  of  Spall  Fractures  in 
2014-T5  Aluminum,  1.6  mm  Impactor 
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Figure  42  Spall  Results  for  2014-T6  Aluminum 
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Figure  44  Results  from  Buffered  Direct  Impact  Shots 

Film  records  from  these  shots,  along  with  two  other  quartz 
gauge  records  discussed  in  the  wave  profile  section  which  fol¬ 
lows,  were  reduced  to  quartz  stress-time  using  "QZ",  a  prograi 

(25) 

for  quartz  gauge  analysis .  Hugoniot  points  were  then  com¬ 

puted  using  appropriate  quartz  stresses  and  impact  velocities 
for  shots  presented  in  Fig  ;res  43  and  44.  The  following  equa¬ 
tions  were  used  for  these  Hugoniot  computations: 
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Results  of  the  above  computations  are  listed  in  Table  V. 

TABLE  V 

HUGONIOT  DATA  FOR  2014-T6  ALUMINUM 


Impact 

2014- 

T6  A1 

Shot 

Velocity 

Quartz  a 

c 

u 

No. 

(mm/usec) 

(kbar) 

(kbar) 

(mm/jjsec) 

9765 

0.117 

9.2 

9.2 

0.056 

9766 

0.311 

24.5 

24.5 

0.149 

9767 

0.261 

9.4 

36.6 

0.2.25 

9768 

0.349 

12.5 

48.7 

0.302 

Probable 

error  is  2%  for  direct  impact 

results  and  3%  for 

buffered 

direct  impact 

results .  This 

is  computed 

as  follows 

Direct  Impact 
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Buffered  Direct  Impact 


*■{$)'  (er*  m’  [(?)••(?) 


—  *  +  l  (1) 2  (.  01) 2  +  (1)2[(.02)2  +  (.01)2]  =  0.025 


VI  \  2/AVl\2 


2Z,  u,. , 
b  Al 


-~^^J  (2)2(.005)2  +  (.25)  2[  (.01) 2  +  (.02) 2  ]  +  (1)2(.01)2  =  0.015 
UA1  " 

Probable  error  =  W  ( . 025) 2  +  ( . 015)  2  J  X  100  =  2.9% 


X  100  =  2.9% 


A  linear  curve  fit  to  hugoniot  data  on  2014-T6  aluminum  re¬ 
sults  in 


a  =  0.37  +  160  u 


for  9  <  a  <  50  kilobars 


(III-27) 


where  o  is  hugoniot  stress  (kilobars) ,  and  u  is  particle  velo¬ 
city  {mm/’^sec)  .  In  Figure  45  this  fit  is  compared  to  measured 
hugoniot  data,  an  acoustic  line  given  by 


a  -  p  CTu  =  179  u 
■o  L 


(III-28) 
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Hugoniot  of  2014-T6  Aluminum  to  50  Kilobars 
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This  hugoniot  may  be  compared  to  the  hugoniot  hydros tat  from 
elastic  constants  (Figure  39) .  For  an  elastic-perfectly  plas¬ 
tic  material,  the  hugoniot  will  lie  ^  1/3  Y  above  the  hydro- 
stat  in  the  stress-particle  velocity  plane,  where  Y  is  the 
yield  stress  measured  under  uniaxial  stress  conditions.  Thus, 
this  amount  may  be  subtracted  from  the  experimental  hugoniot, 
above  the  precursor  elastic  limit,  to  obtain  the  hydrostat. 

Using  appropriate  values  for  2014-T6  aluminum,  we  find  that 
the  experimental  hydrostat  lies  ^  1%  above  the  hydrostat  cal¬ 
culated  from  elastic  constants. 

4.  Wave  Profiles 

Wave  profiles  were  measured  using  both  quartz  transducers  and 
the  laser  velocity  interferometer.  Quartz  data  was  first  re¬ 
duced  to  interface  stress-time  between  the  aluminum  target  and 
the  quartz  gauge.  Interferometer  data  was  first  reduced  to 

interface  velocity-time  using  an  analysis  similar  to  that  re- 

(22) 

ported  by  Barker.  Interface  velocity  and  interface  stress 

were  then  reduced  to  stress  states  in  aluminum  using  standard 
impedance  matching  procedures . ^ 26 ^ 

Transmitted  wave  profiles,  measured  using  quartz,  are  presented 
in  Figure  46  for  a  propagation  distance  of  'v  3  mm  and  impact 
velocities  of  0.125  and  0.432  mm/ysec.  A  ^  1.5  mm  fused  quartz 
impactor  was  used  for  both  of  these  shots.  The  lower  stress 
shot  (9769)  showed  a  peak  stress  level  of  ^  9  kilobars.  This 
peak  stress  is  just  above  the  elastic  limit,  4.5  <  a„_T  <  10  kil- 
obars,  and  no  separation  of  elastic  and  plastic  waves  is  visible 
on  the  record;  however,  there  is  a  substantial  amount  of  ramping 
on  top  of  the  wave.  We  do  not  completely  understand  why  the  rise¬ 
time  is  so  poor  on  shot  (9769),  ^  0.05  ysec.  This  amount  is  about 
twice  the  risetime  that  can  be  accounted  for  due  to  impact  tilt. 
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Transmitted  wave  profiles  measured  using  the  laser  velocity  in¬ 
terferometer  are  shown  in  Figures  47  and  49.  In  Figure  47  pro¬ 
files  are  shown  after  propagation  distances  of  'v  1,  3,  6.5,  and 
12  mm  for  an  impact  stress  of  'v  20  kilobars.  Fused  quartz 
impactors  were  used  for  this  series  which  shows  development  of 
elastic,  plastic,  and  release  waves  in  2014-T6  aluminum  with 
propagation  distance.  Several  significant  features  may  be 
noted:  elastic  precursor  decay  with  propagation  distance, 
spreading  of  the  plastic  wave,  and  spreading  of  the  release 
wa  /e. 


Figure  47  Wave  Profiles  in  2014-T6  Aluminum  for  an 
Average  Impact  Velocity  of  0.284  mm/psec 
and  Propagation  Distances  of  1,  3,  6.5, 
and  12  mm  (Fused  Quartz  Impactors) 


Wave  spreading  with  propagation  distance  is  best  shown  by 
plotting  "phase  lines",  e.g.,  lines  of  constant  stress  in  the 
x-t  plane.  An  x-t  plot.  Figure  48,  was  constructed  from  pro- 
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files  presentd  in  Figure  47  for  both  the  compression  and  rare¬ 
faction  waves.  If  phase  lines  are  straight,  then  wave  propa¬ 
gation  is  steady-state.  Note  that  steady-state  conditions  are 
approached  after  o  12  ran  for  the  compressive  wave  front  as  the 
"phase  lines"  become  parallel,  straight  lines.  Separation  of 
elastic  and  plastic  waves  is  clearly  visible  in  these  loading 
phase  lines;  this  is  shown  as  a  gap  between  phase  lines  in  the 
6  to  10  kilobar  range.  This  gap  broadens  with  propagation 
distance. 

In  Figure  48  the  rarefaction  phase  lines  are  seen  to  fan  out 
from  the  impact  surface.  There  is  slight  evidence  of  an  elas¬ 
tic-plastic  type  release.  This  is  visible  as  a  broadening  of 
phase  line  separation  towards  the  middle  of  this  release  wave 
fan. 

Although  not  carried  out  here,  analytic  forms  of  loading  and 
unloading  phase  lines  are  very  useful  for  predicting  wave 
attenuation  with  propagation  distance.  Shot  9779  in  Figure  47 
is  an  example  of  an  attenuated  wave  profile,  which  occurs  when 
loading  and  release  phase  lines  overlap  in  the  x-t  plane. 

In  Figure  49  the  compressive  wave  is  shown  as  a  function  of 
compressive  stress.  Note  that  increasing  peak  stress  from 
38  to  48  kilobars  changes  the  precursor  very  little.  Further, 
with  this  increase,  the  elastic-plastic  wave  separation  and 
plastic  wave  risetime  are  changed  only  slightly. 

Elastic  precursor  amplitudes  are  summarized  in  Table  VI  for 
2014-T6  aluminum.  Figure  50  is  a  plot  of  these  data  for 
'v  20  kilobar  impact  stress.  Bars  attached  to  these  data  re¬ 
present  the  uncertainty  in  identifying  amplitudes  due  to  the 
rounded  nature  of  the  precursor. 
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Figure  49  Wave  Profiles  in  2014-T6  Aluminum  for  a  Propaga¬ 
tion  Distance  of  'v  3  mm  and  Two  Peak  Stress 
Levels  (Tungsten  Carbide  Impactors) 
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SPALLATION  WAVE  PROFILES 


In  order  to  explain  significant  features  of  measured  wave  pro¬ 
files,  wave  interactions  are  first  illustrated  in  Figure  51. 
Several  simplifying  assumptions  are  made.  The  shock  is  shown 
as  a  single  wave  rather  than  an  elastic-plastic  wave  system, 
and  small  interactions  between  shock  waves  have  been  omitted. 
Fracture  of  the  material  is  assumed  to  occur  along  a  single 
surface  (spall  plane)  rather  than  in  a  volume  around  this  plane 
(Figures  40  and  41) .  Wave  interactions  a-“  the  spall  plane  are 
shown  as  if  the  material  exhibits  a  time  dependent  fracture, 
i.e.,  the  percentage  of  fractured  area  in  the  spall  plane  in¬ 
creases  with  time  until  the  entire  spall  surface  has  separated. 
During  the  early  stages  of  the  fracture  process,  before  complete 
separation  is  achieved,  waves  crossing  the  spall  plane  are  par¬ 
tially  transmitted  and  partially  reflected.  After  complete 
separation  has  occurred,  any  waves  trapped  in  the  spalled  piece 
will  reverberate,  giving  a  decaying,  sinusoidal  motion  to  the 
rear  surface  as  the  wave  is  damped  by  viscous  forces. 
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Figure  51  Wave  Interactions  at  Spall  Plane  and  Resultant 
Velocity-Time  and  Stress-Time  Histories  at 
Free  Rear  Surface 
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Details  of  the  free  surface  motion  of  an  aluminum  impactor 
striking  an  aluminum  target  are  shown  in  Figures  52  and  53. 

The  wave  front  exhibits  elastic-plastic  flow  as  there  is  a 
break  at  about  0.1  mm/ysec  (Figures  52  and  53).  This  corre¬ 
sponds  to  an  elastic  wave  of  approximately  7  kbars  proceeding 
the  plastic  waves.  The  plastic  wave  ramps  towards  a  final 
velocity  which  remains  constant  until  rarefactions  originating 
from  the  free  rear  surface  of  the  impactor  arrive  and  decrease 
the  velocity  of  the  surface.  For  tests  in  which  no  spall  occurs, 
this  decrease  continues  until  the  velocity  is  returned  nearly  to 
zero,  e.g.,  lower  trace  in  Figure  53. 


-  Vj=0.253mm/)Jsec 


Vj=0.194min/ysec. 
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Figure  52  Spall  Wave  Profiles  for  0.6  mm  Impactor  Series 
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TIME— MICROSECONDS 

Figure  53  Spall  Wave  Profiles  for  1.6  mm  Impactor  Series 


Evidence  of  fracture  in  and  about  the  spall  plane  is  shown  as 
a  reversal  in  the  velocity  of  the  surface  at  approximately 
.26  microseconds  in  Figure  52  and  .7  microseconds  on  the  upper 
trace  in  Figure  53.  Signals  from  fractured  surfaces  reach  the 
rear  surface  with  components  of  velocity  opposite  in  sign  to 
the  rarefaction  waves.  As  fracturing  progresses,  the  surface 
is  no  longer  "pulled  back"  by  the  release  wave  system,  but  ac¬ 
celerates  ,  only  to  be  decelerated  again  as  the  entrapped  wave 
reverberates  back  and  forth  between  the  spall  surface  and  the 
specimen  rear  surface.  If  it  is  assumed  that  the  decrease  of 
the  free  surface  velocity  to  the  point  of  reversal  (referred 
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to  as  "pullback")  is  related  to  the  maximum  tensile  stress  at 
the  spall  plane,  a  quantitative  measure  of  the  negative  stress 
required  to  fracture  the  material  for  a  given  wave  shape  is  pro¬ 
vided  by  the  profile.  The  functional  form  of  the  relationship 
between  the  maximum  tensile  stress  and  pullback  has  not  been 
firmly  established,  but  attempts  to  correlate  experimental  re¬ 
sults  with  calculations  of  spall  plane  stress  have  resulted  in 
the  following: 

ot  =  pC  §£  (IH-2S) 

where  a  is  maximum  tensile  stress  at  the  spall  plane,  p  and 
C  are  local  density  and  dilatational  wave  velocity,  and  AV  is 
pullback . 

Using  this  relation  gives  a  'v  12  kilobars  for  the  thin  impactor 
series  and  'v  15  kilobars  for  the  thick  impactor  series.  This 
behavior  is  backwards  from  results  reported  earlier  for  6061-T6 
aluminum.  For  6061-T6  aluminum  pullback  was  smallest  for  the 
thinner  impactor  series.  Reasons  for  this  difference  are  not 
understood  at  this  time. 

COMPARISON  OF  SHOCK  PROPAGATION  AND  FRACTURE 
CHARACTERISTICS  OF  2014-T6  AND  6061-T6  ALUMINUM 

1.  Elastic  Constants 

The  elastic  constants  and  derivatives  of  elastic  constants  of 
2014-T6  aluminum  and  6061-T6  aluminum  are  compared  in  Table  VII. 
In  Figure  54  comparison  is  made  between  the  hugoniot  hydrostat 
for  these  two  alloys  as  calculated  from  Eqs .  III-19,  III-20, 
III-21  and  III-22  using  some  of  the  values  listed  in  Table  VII. 
These  calculated  hugoniots  differ  by  v  3%  at  50  kilobars. 
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Figure  54  Hugoniot  Hydrostats  Calculated  for  2014-T6  and 
6061-T6  Aluminum  Using  Eqs .  (III-19) ,  (III-20)  , 
(III-21) ,  and  (III-22) (29) 
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2.  Spall  Threshold 

Within  experimental  uncertainty,  there  is  no  difference  be¬ 
tween  spallation  data  obtained  by  recovery  tests  on  2014-T6 
aluminum  and  published  results  on  6061-T6  aluminum  (Figure  55) . 
Spallation  data  available  on  6061-T6  aluminum  might  be  used 
with  good  accuracy  to  predict  the  spallation  behavior  of 
2014-T6  aluminum. 


3.  Hugoniot 

In  Figure  56  experimental  hugoniots  for  2014  and  6061  are  com¬ 
pared.  The  hugoniot  for  2014  lies  above  the  hugoniot  for 
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6061  by  about  4%  at  50  kilobars.  There  is  a  3%  difference  be¬ 
tween  calculated  hugoniots  at  50  kilobars  (Figure  54) .  This 
suggests  that  part  of  the  difference  between  experimental  hug¬ 
oniots  is  due  to  differences  in  the  elastic  constants  and  den¬ 
sity  of  the  two  materials.  Any  difference  in  precursor  ampli¬ 
tude  would  also  contribute  to  the  difference. 


Figure  56  Hugoniot  for  2014-T6  Aluminum  Compared 
to  Hugoniot  for  6061-T6  Aluminum(29) 
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4.  Elastic  Precursor  Decay 

Elastic  precursor  decay  observed  for  2014-T6  aluminum  is  com¬ 
pared  to  6061-T6  aluminum  in  Table  VIII.  Values  listed  for 
6061-T6  aluminum  represent  averages  over  various  plate  stocks 
that  have  been  studied  by  General  Motors  Corporation. 
Precursor  amplitudes  and  decay  rates  agree  within  the  uncer¬ 
tainty  of  defining  o„_T . 


TABLE  VIII 


COMPARISON  OF  PRECURSOR  AMPLITUDE  FOR  2014-T6  AND  6061-T6 
ALUMINUM  AT  ^  20  KILOBARS  IMPACT  STRESS  ' 


Propagation 

™  i  i  i  j  . . 

Precursor 

Amplitude 

Distance 

2014-T6  Al* 

6061-T6  Al** 

(mm) 

(kbar) 

(kbar) 

'v  3 

6.5  ±  .5 

6.5 

^  6.5 

5.5  ±  .5 

5.5 

%  12 

4.5  ±  .5 

4.5 

*  12.5  mm  plate  stock. 


*  Average  value  from  3.2,  6.4  and 
12.5  mm  plate  stocks. 
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SECTION  IV 

CONCLUSIONS 


SUBSEQUENT  YIELD  STUDY 

1.  The  yield  surfaces  are  convex. 

2.  Cross-effects  do  exist. 

3.  There  is  evidence  of  a  substantial  Bauschinger  effect. 

4.  The  subsequent  yield  surface  both  expands  and  translates, 
indicating  that  the  hardening  mechanisms  are  more  compli¬ 
cated  than  both  the  isotropic  and  kinematic  hardening. 

5.  The  general  hardening  theory,  attributed  to  Hodge,  best 
describes  the  subsequent  yield  behavior  of  the  material. 

6.  No  evidence  of  corners  was  found  in  the  study;  however 
from  the  definition  of  yield  used,  corners  would  not  be 
expected. 


k 


HEATING  RATE  STUDY 

1.  At  500°F,  yield  and  ultimate  strength  is  most  sensitive 
to  heating  rates  between  1.0°F/sec  and  0 . 01°F/second. 

The  heating  rate  itself  does  not  appear  to  be  important? 
the  time  above  the  artificial  aging  temperature  is  most 
significant. 
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2.  The  strength  at  500°F  when  heated  to  temperature  in  less 
than  30  minutes  (0.25°F/sec)  is  approximately  three  times 
the  strength  as  when  50  hours  is  required  to  reach  tempera¬ 
ture  (0 . 0025°F/second) . 

3.  Above  a  heating  rate  of  0.25'’F/sec  heating  rate  effects, 
if  any,  are  lost  in  the  data  scatter  and  the  predominant 
influence  on  strength  is  the  test  temperature. 

4.  Coalescence  of  precipitates  appears  to  be  the  major  mech¬ 
anism  influencing  the  change  of  material  properties  with 
heating  rate.  Therefore,  different  precipitation  rates 
would  influence  the  effect  of  heating  rate  on  aluminum 
alloys . 

5.  Little  effect  on  precipitation  is  observed  metallographi- 
cally  at  heating  rates  above  0.25°F/sec  (30  minutes  to 
reach  500°F) ,  and  the  alloy  strength  appears  to  be  insen¬ 
sitive  to  time  except  for  the  loss  due  to  temperature 
itself  ('v  30%  for  yield  at  500°F)  . 

SHOCK  WAVE  PROPAGATION  STUDY 

1.  The  shock  wave  response  of  2014-T6  aluminum  is  elastic- 
plastic  and  nearly  identical  to  that  of  6061-T6  aluminum 
except  for  a  difference  between  experimental  hugoniots 

4%  at  50  kilobars) . 

2.  The  elastic  precursor  amplitude  ranged  from  8  kilobars 
after  a  propagation  distance  of  1  mm  to  4.5  kilobars 
after  12  mm. 
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3.  The  hugoniot  for  2014-T6  aluminum  may  he  approximated  by 

Stress  Range  Accuracy 

a  =  179  u  0<a<8  %  1% 

a  =  0.37  +  160  u  8  <  a  <  50  %  3% 

where  a  is  stress  in  kilobars  and  u  is  particle  velocity 
in  mra/ysec. 

4.  For  the  Murnaghan  equation  of  state,  the  hugoniot  calculated 
from  elastic  constants  lies  above  experimental  mean  stress 
hugoniot  by  v  3%  at  50  kilobars. 

5.  The  spall  threshold  for  2014-T6  aluminum  is  dependent  upon 
loading  time.  Impact  velocity  (aluminum  impactor)  necessary 
for  incipient  spall  decreased  from  0.216  mm/ysec  for  0.19 
ysec  pulse  duration  to  0.180  mm/ysec  for  0.5  ysec  pulse  dur¬ 
ation. 

6.  For  a  peak  stress  of  20  kilobars  the  compressive  wave  front 
becomes  steady  after  a  propagation  distance  of  %  12  mm. 

For  this  same  peak  stress  the  release  wave  is  a  simple 
centered  rarefaction  fan  with  only  slight  evidence  of 
elastic-plastic  structure. 
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APPENDIX  A 

MATERIAL  DESCRIPTION 


Three  different  aluminum  alloys  were  used  in  this  program; 
6061-T6  and  2014-T6  aluminum  alloys  were  used  for  the  high 
heating  rate  studies  (Section  II) ,  the  same  2014-T6  alloy 
for  the  shock  wave  study  (Section  III) ,  and  2024-T3  aluminum 
alloy  for  the  biaxial  studies  (Section  I) . 

The  6061-T6  material  was  obtained  in  the  form  of  0.25  inch 
diameter  rod.  Photomicrographs  of  the  material  (Figure  A-l) 
show  elongated  grains  running  parallel  to  the  rod  axis.  The 
2024-T3  alloy  obtained  in  the  form  of  3.0  inch  diameter  rod, 
shews  a  similar  structure  but  with  grains  more  severely  elong¬ 
ated  (Figure  A-2) .  Both  of  these  materials  exhibit  little 
grain  size  variation  across  the  diameter  of  the  rod. 

The  2014-T6,  obtained  in  the  form  of  0.5  inch  thick  plate, 
exhibited  non-uniform  grain  size  through  the  thickness  of  the 
plate.  The  grain  size  decreased  from  the  plate  outer  eage  to 
the  plate  center  (see  Figure  A-3) . 

Certified  chemical  composition  and  heat  treatment  of  the 
materials  are  listed  in  Table  A-I. 
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TABLE  A- I 

CHEMICAL  COMPOSITION 


6061-T6 


MIN. 

MAX. 

(%) 

(%) 

Si 

0.4 

0.8 

Fe 

0.7 

Cu 

0.15 

0.40 

Mn 

0.15 

Mg 

0.8 

1.2 

Zn 

0.25 

Cr 

0.15 

0.35 

Ti 

0.15 

Other  (Each) 

0.05 

Others  (Total) 

0.15 

AL 

REMAINDER 

Heat  Treatment 

T6- 

51 

2014-T6 

2024- 

-T3 

(NOMINAL) 

MIN. 

MAX. 

(%) 

(%) 

(%) 

0.8 

0.50 

0.50 

4.4 

4.9 

0.8 

0.3  0.9 

0.5 

1.2  1.8 

0.25 

0.1 

0.05 

0.15 

remainder 

REMAINDER 

T6-51 

T3-51 
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